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ABSTRACT 


A  theoretical  and  experimental  investigation  of  thrust 
control  by  secondary  gas  injection  has  been  completed. 
Thrust  control  is  treated  from  the  two  aspects  of  thrust 
modulation  and  vector  control  with  the  emphasis  being  on 
the  former. 

The  thrust  modu^latlonT study  is  concerned  with  the 
use  of  injection  as  a  means  of  throttling  the  main  nozzle 
flow.  Since  experimental  data  shows  that  thrust  is  propor¬ 
tional  to  total  nozzle  flow,  thethrustis  thus  also  proportional 
to  flow  throttling  or  injection  rate.  The  performance  (flow 
throttling  and  thrust)  of  avariety  of  secondary  injection  con¬ 
figurations  was  measured  with  a  low  pressure-ratio,  axially 
symi-netric  nozzle  as  well  as  a  high  pressure-ratio,  two- 
dimensional  nozzle.  Cold  air  (530 *R)  was  the  working  fluid 
in  the  primary  and  secondary  flows.  The  flow  throttling 
results  agree  with  an  analysis  based  on  three  flow  models 
and  the  schlieren  photographs  (from  the  two-dimensional 
tests)  show  that  all  three  types  of  flow  occur.  The  effects 
of  the  important  variables  are  determined  analytically 
and/or  experimentally. 

The  vector  control  study  concerns  asymmetric 
secondary  gas  injection  as  a  means  of  controlling  the  direc¬ 
tion  of  the  thrust  vector.  Some  experimental  data  was  ob¬ 
tained  from  the  two-dimensional  nozzle  tests.  This  data 
is  compared  to  data  from  other  sources  and  an  explanation 
is  given  for  the  process  which  governs  the  generation  of 
side  force. 
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SECTION  I 
INTRODUCTION 

The  wide  variety  of  missions  presently  associated  with  vehicles 
using  large  rocket  engines  has  made  apparent  the  need  for  versatility 
in  engine  performance.  Of  particular  importance  in  this  respect  is  the 
capability  of  rocket  thrust  control.  Although  much  work  has  been  done 
on  this  subject,  the  state  of  the  art  certainly  does  not  preclude  the  pos¬ 
sibility  of  further  improvements.  The  conventional  gimbal  system, 
especially  in  the  largest  engines,  is  enormously  complicated  and  not  the 
ultimate  in  efficiency.  In  the  case  of  solid  propellant  engines,  the  prob¬ 
lem  of  thrust  control  is  particularly  critical  in  view  of  the  inherently 
"uncontrollable"  nature  of  the  burning  process.  Liquid  propellant  devices, 
however,  are  not  as  rigidly  fixed  in  their  basic  mode  of  operation  and  thus 
may  be  better  suited  to  modifications  designed  to  incorporate  thrust  con¬ 
trol  capability.  In  either  case,  the  need  for  more  efficient  systems  is 
great. 

The  Naval  Supersonic  Laboratory  (NSL)  has  conducted  an  investi¬ 
gation  of  a  promising  method  of  thrust  control.  The  method  is  one  in  which 
a  secondary  flow  is  injected  into  a  main  nozzle  to  change  the  magnitude 
and/or  direction  of  the  nozzle  flow.  Since  nozzle  thrust  is  directly  pro¬ 
portional  to  the  product  of  mass  flow  and  velocity,  any  change  in  these 
quantities  implies  a  change  in  thrust.  The  present  study  deals  with  sym¬ 
metric  injection  for  flow  throttling  (thrust  modulation)  and  asymmetric 
injection  for  flow  direction  control  (vector  control).  Most  of  the  work, 
however,  concerns  symmetric  inje-etion  configurations. 

A  variety  of  parameters  have  been  investigated  in  the  two  phases 
of  research  which  constitute  the  overall  program.  The  Phase  I  portion 
centered  around  three-dimensional  air  flow  tests  and  included  analytical 
studies.  This  work  is  reported  in  detail  in  Ref.  1  and  summarized  briefly 
herein.  Two-diniensional  air  flow  tests  formed  the  basis  of  the  Phase  II 
effort  which  is  described  in  detail  in  this  report. 


TR  448 


1 


CONFIDENTIAL 


CONFIDENTIAL 


SECTION  II 
THRUST  MODULATION 

Analytical  and  experimental  studies  have  been  carried  out  which 
investigated  symmetric  secondary  injection  as  a  means  of  thrust  modu¬ 
lation,  The  first  portion  of  this  section  briefly  summarizes  the  most 
important  features  of  the  work  accomplished  in  the  Phase  I  program. 

The  Phase  I  program  inrliiderl  hoth  analytical  and  e-xperimental  studies- 
The  second  portion  of  this  section  presents  further  analytical  work.  This 
involves  an  analogy  between  the  secondary  injection  and  blunt  body  problems 
as  well  as  extensions  of  the  Phase  I  analytical  developments.  A  presen¬ 
tation  of  the  details  of  the  Phase  II  two-dimensional  air  flow  tests  con¬ 
stitutes  the  third  part  of  this  section  while  the  fourth  and  final  one  is  a 
correlation  and  comparison  of  NSL  and  other  data. 

A.  PHASE  I  -  THEORY 

The  first  thorough  attempt  to  find  an  analytical  solution  to  the 
flow  throttling  problem  was  made  by  A.  I.  Martin  in  Ref.  2.  Martin  re¬ 
garded  secondary  injection  as  a  possible  method  of  obtaining  a  variable- 
throat  nozzle  for  use  in  controlling  the  performance  of  a  jet  engine  at 
different  operating  conditions.  Two  theoretical  models  were  developed 
although  this  procedure  was  adopted  merely  as  a  simplification  necessary 
to  generate  an  amenable  analytical  problem.  Martin  expected  tb.at  the 
real  case  would  be  a  combination  of  the  features  of  both  models. 

The  problem  solutions  obtained  by  Martin  came  from  appropriate 
applicati'^'M  of  simplified  forms  of  the  equations  of  fluid  motion.  He  com¬ 
pared  numerical  results  to  some  test  data  and  the  agreement  between 
theory  and  experiment  was  not  very  good.  In  addition,  appropriate  use 
of  Martin's  solutions  gave  results  which,  like  his  own  data,  disagreed 
with  NSL  experimentai  data  (reported  in  Section  B). 

Further  study  of  the  Martin  analyses  led  to  the  development  of 
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several  NSL  analytical  models.  There  are  three  NSL  flow  models  of 
which  two  are  modifications  of  the  Martin  ''es.  These  are  the 
"secondary  mixing"  (modified  "mixlnf;  process")  model  and  the  "sheet 
flow"  (modified  "vortex  sheet")  nior-  /.n  entirely  new  one  has  also 
been  developed  called  the  "secondary  expansion"  model.  The  most  im¬ 
portant  features  of  the  NSL  analyses  are  described  below  while  a.  de¬ 
tailed  treatment  may  be  found  in  Ref.  1. 

Although  the  theoretical  flo-.v  configurations  are  treated  separately 
below,  some  simplifying  assumptions  are  common  to  all  of  them.  The 
most  basic  one  is  that  the  flow  is  one -dimensional.  In  addition,  the 
primary  and  secondary  fluids  are  assumed  to  be  perfect  gases  whose  flow 
is  isentropic.  These  assumptions  may  not  be  satisfied  but  their  use  leads 
to  an  approximately  correct  choice  of  dependence  of  problem  variables. 

The  assumption  of  one -dimensional,  perfect  gas  flow  immediately 
relates  certain  problem  variables.  The  equation  of  state  for  a  perfect 
gas, 


p  =  p  R  T 


(1) 


relates  the  pressure,  density,  and  temperature  at  any  point  in  the  flow. 
As  shown  in  Ref.  3,  the  energy  equation  may  be  written  in  the  following 
appropriate  form: 


J 

To 


(2) 


The  energy  equation  is  also  useful  in  obtaining  the  following  one¬ 
dimensional  isentropic  flow  relations: 


TR  448  4 

CONFIDENTIAL 


CONFIDENTIAL 


All  of  the  above  equations  apply  throughout  the  flow  in  all  of  the 
theoretical  models. 

1  Secondary  Mixing  Model 

In  the  secondary  mixing  case  the  secondary  flow  mixes  with  the 
primary  stream  to  produce  a  combined  stream.  The  important  features 
of  the  flow  configuration  are  shown  in  the  sketch  below. 


The  secondary  flow  (rate  w^)  enters  the  nozzle  (radius  Rj  )  through  a 
sonic  orifice  (area  and  is  directed  at  an  angle  (})  with  respect  to 

the  nozzle  axis.  The  primary  and  secondary  streams  mix  to  form  a 
final  stream  tube  of  radius  .  The  primary  stream  (rate  w^)  flows 
uniformly  at  Mach  number  Mj  at  station  1  while  at  station  2  a  uniform 
combint^d  stream  flows  at  Mach  number  . 

The  equations  of  continuity  and  momentum  are  applied  to  the 
flow  and,  in  particular,  are  evaluated  at  stations  1  and  2.  The  continuity 

TR  448  5 


CONFIDENTIAL 


CONFIDENTIAL 


equation  may  be  written  as 


r  1 

r  1 

2 

r  1 

£2. 

5i 

M> 

li 

£i 

M, 

II 

Po 

Ri 

To 

Po 

To 

■p-^s' 

POs’ 

[T’al 

Po 

"TT 

(5) 


where  the  primary  and  secondary  fluids  are  the  same  gas. 
The  momentum  equation  is  given  by 


£j. 

Po 


T 

.  *  .  i 

F>i>8 

'P»s  ' 

1  +  Y  Ml 

-  V 

It  R,' 

Po 

cos 


£2 

Po 


1  +  Y  Mi' 


(6) 


in  which  the  primary  and  secondary  gases  are  again  the  same.  This 
equation  has  the  additional  assumption  that  the  stagnation  conditions  of 
the  combined  stream  are  equal  to  those  of  the  primary  stream.  Equa¬ 
tion  (6)  also  has  a.",  assumption  which  concerns  the  "stagnant  air"  zone 
shown  in  the  sketch  on  page  5.  It  was  felt  that  this  separated  flow  zone 
might  exist  at  some  very  low  pressure.  That  is,  the  centrifugal  force 
on  the  ;jecondary  flow  as  it  changes  direction  by  (it  -  (^)  might  leave 
little  flow  near  the  wall  just  downstream  of  the  injection  point.  There¬ 
fore,  for  the  purposes  of  analysis,  it  was  assumed  that  the  pressure  in 
this  zone  is  negligible  (mathematically,  zero)  compared  to  that  in  the 
primary  and  secondary  streams. 


All  of  the  work  in  this  report  concerns  configurations  where  £  90*. 
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It  should  be  pointed  out  that  Eqs.  (5)  and  (6)  are  not  in  a  form 
■which  corresponds  to  the  general  representation  of  the  problem.  This 
is  due  to  the  assumption  that  the  primary  and  secondary  flews  are 
the  same  gas.  The  case  of  different  gases,  however,  does  not  permit 
the  use  of  the  sirr.piified  equations  shown  above.  The  values  of  the 
gas  constant  R  and  the  ratio  of  specific  heats  y  must  be  determined  for 
the  combined  stream  from  trial  values  of  the  primary  flow  associated 
with  a  given  secondary  flow.  This  iterative  procedure  is  thus  super¬ 
imposed  on  the  remaindv;r  of  the  method  of  solution.  Apart  from  this 
matter,  the  method  of  solution  is  complicated  above  and  beyond  that 
of  the  ”  same  gas"  case  due  to  another  factor.  This  factor  is  that  the 
mixing  of  gases  implies  an  increase  in  entropy  (which  also  means  a 
loss  in  efficiency).  Thus  Eqs.  (3)  and  (4)  (and  thus  (5)  and  (6))  would 
be  replaced  by  the  more  complicated  relationships  appropriate  to 
the  situation.  Hc.wever,  here  again  the  procedure  is  an  iterative 
one  since  the  entropy  increase  is  a  function  of  the  amount  of  each  gas 
or,  in  effect,  the  primary  flow  associated  with  a  given  secondary  flow. 

It  is  fairly  obvious  at  this  point  that  the  "different  gas"  case  thus  be¬ 
comes  indeed  complicated.  In  fact,  its  cc'mplexity  is  out  of  proportion 
to  that  of  the  other  models  as  well  as  its  usefulness  as  an  approximate 
analysis.  Furthermore,  by  way  of  experience  with  data,  it  is  difficult 
to  justify  in  view  of  the  fact  that  the  "same  gas"  case  is  the  poorest 
(compared  to  data)  of  the  theoretical  models.  Therefore,  the  remainder 
of  this  discussion  concerns  the  "same  gas"  case  and  all  references  to 
the  secondary  mixing  model  pertain  to  this  case  only. 

The  flow  configuration  sketch  on  page  5  shows  that  th(?  secondary 
flow  is  injected  through  a  .sonic  orifice  of  area  A,^^  from  stagnation  ; 
pressure  and  temperature  po^  and  ,  respectively.  In  the  form 
of  a  non-dimensional  ratio  which  compares  it  to  the  original  (wg  =  0) 
primary  flow,  the  secondary  flow  rate  is 
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where  the  two  flows  are  the  same  gas.  Equation  (7)  is  the  special 
case  (primary  gas  =  secondary  gas)  of  the  general  form  given  in  the 
iiext  section  (A2).  The  choice  of  this  parameter  (ws/wp^)  as  an  index 
of  the  secondary  flow  rate  is  somewhat  arbitrary  although  its  use¬ 
fulness  is  made  apparent  by  the  discussion  in  the  next  paragraph. 

When  certain  final  assumptions  are  made  the  continuity  and 
momentum  equations  ((5)  and  (6))  yield  a  solution  to  the  problem.  The 
geometry  is  specified  by  the  assumption  of  values  for  the  parameters 
A^^/ttRj*  and  ^  while  the  particular  gas  employed  designates  The 
final  assumption  necessary  is  an  overall  pressure  ratio  (pj/po)  which 
implies  the  value  of  the  Mach  number*  Mj  .  Since  the  pressure, 
density,  and  temperature  ratios  are  functions  of  Mach  number 
(according  to  Eqs.  (2),  (3),  and  (4)),  Eq8(5)  and  (6)  thus  become  a 
pair  of  simultaneous  equations  with  unknowns  [R2/Ki]^  and  Mj  . 

The  particular  procedure  used  in  solving  the  equations  is  described 
in  Appendix  A.  As  shown  in  this  appendix,  flow  throttling  results 
are  conveniently  handled  in  terms  of  non-dimensional  ratios.  Thus, 
the  results  for  a  typical  case  are  shown  in  Fig.  1  ,  taken  from  Ref.  1 , 
where  primary  flow  ratio  is  plotted  against  secondary  flow  ratio.  For 
purposes  of  comparison,  also  shown  is  the  Martin  mixing  process 
result. 

It  should  be  noted  that  in  Fig.  1  as  well  as  subsequent  figures 
the  flow  throttling  curves  do  not  begin  at  Ws/wp^  =  0.  This  results 
from  the  assumption  of  sonic  injection  of  secondary  flow.  The  lower 
limit  of  wg/wjjg  is  thus  determined  according  to  the  assumption  that 
the  sonic  secondary  pressure  is  approximately  the  same  as  the  sonic 
primary  pressure.  In  equation  form. 


Po 


Pi 


*  Results  show  that  maximum  flow  throttling  occurs  when 

^  (Pz/Pu  =  0-  5283).  This  is  understandable  since  the  primary 
stream  occupies  a  minimum  of  area  at  sonic  conditions.  In  all  of 
the  theoretical  calculations,  therefore,  Mj  is  unity. 
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which  determines  the  minimum  POg/Po  •  Since  Pz/Po  is  the  sonic 
(Mj  =  1)  pressure  ratio  for  air,  if  the  secondary  gas  is  also  air, 
then  the  minimum  Po^/Po  is  unity. 

2.  Sheet  Flow  Model 

The  basic  assumption  in  the  sheet  flow  model  is  that  the 
primary  and  secondary  flows  do  not  mix.  The  flow  configuration  is 
sketched  below. 


^0*  ’^<>8 


\  MM  aMMB  MB 

Stognonf  Air  I 


I  Main  Sifreom 
I  Cylindrical  Channel 


The  secondary  flow  enters  the  nozzle  in  the  same  manner  as  that  in 
the  secondary  mixing  model-  In  the  present  model,  however,  the 
secondary  stream  flows  in  a  sheet  circumjacent  to  the  primary  one 
and  there  are  two  distinct  streams  at  station  2.  At  this  station,  the 
secondary  stream  occupies  an  area  while  it  flows  at  Mach  number 

M^^.  The  pressure  in  the  stagn.'nt  air  zone  is  again  assumed  to  be 
negligible  (mathematically,  ?ero)  compared  to  that  in  tlie  primary 
and  secondary  str.?ams. 
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The  equations  of  continuity  and  momentum  are  applied  to  the 
flow  ati  in  the  secondary  mixing  model.  Now  there  are  two  continuity 
equations  written  ns 


Po 


p  -l 

l 

r  ^ 

1  1  1 

r  1 

Ra 

Ri 

Ml 

To 

1 

II 

2 

II 

To 

(8) 


for  the  primary  flow,  and 
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for  the  secondary  flow. 

The  momentum  equation,  after  substitution  of  Eq.  (9)|  is 
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where  the  primary  and  secondary  streams  may  be  differejit  gases.  T 
area  ratio  isentropic  area  ratio  which  corresponds  to 

where  is  implied  by  the  assumption  that  the  static  pressures 

in  both  streams  are  one  and  the  same  iPz^  ~  Pz)-  This  is  explained 
in  Appendix  B. 

Tile  general  form  of  Eq.  (7)  for  the  secondary  flow  ratio  is 


(11) 
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where  the  primary  and  secondary  gases  may  be  different. 

The  final  assumption  of  values  for  A^^/itRi*  and  ^  as  well  as 
the  gases  involved  leads,  as  in  the  secondary  mixing  model,  to  a 
final  set  of  simultaneous  equations  in  the  unknowns  [_Rj  /R|  ]A  and  . 
Appendix  A  describes  the  graphical  technique  used  to  solve  these 
equations.  Typical  results  are  shown  in  Fig.  1  (taken  from  Ref.  1) 
where  the  comparable  Martin  vortex  sheet  model  has  been  added. 

3.  Secondary  Expansion  Model 

In  the  secondary  expansion  model  the  flow  behavior  is  somewhat 
different  from  either  of  the  preceding  cases.  The  secondary  flow  enters 
the  nozzle  in  the  same  way  and,  as  in  the  sheet  flow  model,  does  not 
mix  with  the  primary  stream.  However,  in  the  present  model  the 
secondary  stream  expands  as  it  flows  downstream,  as  shown  in  the 
sketch  below,  filling  the  region  between  the  wall  and  the  primary  stream. 


This  secondary  expansion  is  isentropic  and  the  Mach  number  Mj  is 
given  by  the  area  ratio*  A^^/A^^^  rather  than  being  implied  by  the 
assumption  that  pj  ~  secondary  expansion  model  assumes 

*  See  the  first  part  of  Appendix  B, 
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that  the  flow  does  not  separate  from  the  wall. 

The  continuity  equations  for  the  present  model  are  Eqs.  (8)  aind  (9) 
while  the  momentum  equation  is 


where  trial  values  of  (which  determines  are  used  according 

to  the  iterative  procedure  explained  in  Appendix  C.  Equation  (12)  takes 
account  of  the  fact  that  the  primai*y  and  secondary  gases  may  not  be 
the  same  as  does  £q.  (I  I)  which  is  the  secondary  flow  ratio  appropriate 
to  the  present  model. 

The  final  specification  of  problem  geometry  and  primary  and 
secondary  gases  again  yields  the  simultaneous  equations  (in  unknowns 
[Rj/Rj]*  and  Mj)  which  were  solved  by  the  graphical  technique  in 
Appendix  A.  The  results  for  a  typical  case  are  again  shown  in  Fig.  1 
(taken  from  Ref.  1). 

B.  PHASE  I- -EXPERIMENT 

A  series  of  experiments  were  performed  as  part  of  the  Phase  I 
program  in  which  secondary  injection  was  used  in  an  axially  symmetric 
nozzle.  A  variety  of  injection  configurations  were  tested  with  air  as 
the  working  fluid  in  both  the  primary  and  secondary  streams.  The 
objective  of  these  tests  was  to  determine  some  of  the  flow  throttling 
and  attendant  thrust  control  characteristics  peculiar  to  the  injection 
method.  A  detailed  treatment  of  this  work  io  reported  in  Ref.  1  while 
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a  brief  summary  Is  presented  in  this  section. 

Air  for  the  primary  and  secondary  flows  was  obtained  from 
au3dliary  wind  tunnel  equipment  at  the  NSL.  F4ow  rate,  stagnatHja  - 
pressure,  and  stagnation  temperature  were  measured  in  both  flows. 

The  primary  stagnation  pressure  was  100  psia  in  nearly  all  of  the 
tests  and  held  constant  regardless  of  the  amount  of  primary  flow  throttling 
Model  forces  were  measured  directly  by  means  of  strain  gage  balance 
systems . 

The  axially  symmetric  nozzle  used  in  these  tests  was  conical 
having  a  subsonic  cone  semi-vertex  angle  of  30°  and  a  supersonic  cone 
semi-vertcxangle  of  15°.  The  design  flow  rate  was  one  pound  per 
second  with  a  stagnation  pressure  of  100  psia  and  stagnation  temperature 
of  530°R.  The  exit  Mach  number  of  the  nozzle,  according  to  the  exit- 
to-throat  area  ratio,  was  2.  1.  Secondary  flow  rates  varied  up  to  about 
0.  15  pounds  per  second  corresponding  to  a  secondary  stagnation  pressure 
of  about  300  psia. 

Secondary  flow  was  injected  from  various  orifice  configurations. 
These  were  usually  holes  equally  spaced  around  the  circumference  of 
the  nazzi  j  throat.  In  one  case  a  slot  replaced  the  holes.  Injection 
angKs  (<|))  of  60°  and  90°  were  tested.  The  circumferential  injection 
angle  P  was  also  varied. 

Some  typical  flow  throttling  data,  taken  from  Ref.  1,  is  presented 
in  Fig.  2  where  it  is  compared  to  the  results  for  the  sheet  flow  and 
secondary  expansion  models.  Configuration  F  has  32  holes  0.032  inches 
in  diameter  while  configuration  J  has  a  slot  0.  020  inches  in  width.  The 
analytical  calculations  are  based  on  the  assumption  that  the  flow  throttling 
takes  place  at  or  near  the  nozzle  throat.  Thus  the  throat  area  (0.430  in*) 
is  assumed  to  be  ttRj*  .  In  addition,  since  the  analytical  methods  do  not 
contain  an  orifice  coefficient*  (C)  in  the  secondary  mass  flow  equation, 
the  area  ratio  was  replaced,  in  the  analytical  calculations,  by 

CA,^^/itRi*.  This  procedure  was  adopted,  somewhat  arbitrarily,  so  that 
the  theoretical  and  experimental  stagnation  pressure  ratios  (Po^/po) 
would  be  about  the  same  for  a  given  wg/wp^. 

*  An  approximation  to  the  average  value  of  C  was  used  for  each 
configuration. 
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The  agreement  between  theory  and  experiment  shown  in  Fig.  2 
i.-?  typical  of  all  of  the  Phase  I  data'^.  These  models  encompassed  almost 
all  of  the  data  and  this  pattern  seems  reasonable  since,  in  a  sense,  they 
are  "extremes"  of  a  type  of  flow  behavior-  In  the  sheet  flow  model  the 
secondary  pressure  (at  station  2  in  the  analysis)  is  determined  by  the 
primary  stream  and,  since  the  change  in  flow  direction  (120°)  is  large, 
the  secondary  stream  is  separated  from  the  nozzle  wall.  In  the  secondary 
expansion  model,  however,  the  secondary  stream  expands  to  fill  the 
area  between  the  primary  stream  and  the  wall.  In  this  case  the  secondary 
pressure  is  independent  of  the  adjacent  primary  value  and  there  is  no 
separated  flow  region.  Some  crude  Phase  I  water  table  observations 
indicated  that  both  flow  configurations  exist  and,  as  the  air  data  seems 
to  prove,  this  may  actually  be  the  case. 

Thrust  control  is  directly  proportional  to  total  noiizle  flow  with 
symmetric  or  asymmetric  injection.  In  over  97  percent  of  the  test 
runs  the  nozzle  thrust  ratio  (F/Fo)  divided  by  the  total  flow  ratio 
[  (Wp  +  Wjj  varies  between  0.96  and  1.04.  This  indicates  that 

the  exhaust  velocity  does  not  cis.ngc  appreciably  with  injection  rate 
and  the  flow  may  behave  as  tltough  the  throttling  had  been  accomplished 
by  an  "effective  geometric  throat". 

C.  FURTHER  ANALYTIC AL  STUD Y- -FLOW  FIELD  ANALOGY 


The  NSL  theoretical  models  (sheet  flow  and  secondary  expansion) 
seem  to  correlate  the  data  fairly  well.  How'ever,  due  to  the  extrinsic 
nature  of  such  analyses  (consideration  of  the  flow  through  a  control 
volume)  only  a  vague  picture  is  acquired  of  the  actual  flow  field.  A 
ne.xt  step  would  seem  to  be  the  development  of  an  analysis  which  gives 
a  clear  understanding  of  the  detailed  nature  of  the  flow  process.  An 
effort  has  been  made  to  formulate  such  an  analysis  of  the  flow  field 
induced  by  secondary  injec*i.on  into  a  supersonic  stream. 

The  water  table  flow  observations,  mentioned  in  connection  with 
tile  theoretical  models,  showed  that  the  main  nozzle  flow  with  injection 
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resembled  the  flow  past  a  blunt  body.  This  suggested  that  if  the  primary 
stream  does  not  mix  greatly  with  the  secondary  stream  perhaps  the 
primary  flow  is  analogous  to  a  blunt  body  flow. 

The  "flow  field  analogy"  analysis  is  formulated  in  two  main 
parts.  First,  the  primary  flow  field  is  treated  in  accordance  with  the 
approximate  method  of  L.i  and  Geiger  (Ref.  4)  for  finding  the  flow  field 
in  the  stagnation  region  of  a  blunt  body.  Second,  the  flow  field  of  the 
secondary  stream  is  formulated  in  accordance  with  similar  consider¬ 
ations.  The  two  flow  fields  are  then  combined  according  to  certain 
conditions  which  must  hold  at  the  boundary  which  separates  the  streams. 
Since  the  combined-flow  problem  is  non-linear,  a  solution  must  be 
obtained  by  approximate  numerical  methods.  Although  a  solution  was 
not  obtained  the  work  accomplished  is  nevertheless  felt  to  be  worthy 
of  mention. 

1.  Li  and  Gieger  Blunt  Body  Analysis 

Since  the  analysis  of  the  primary  flow  field  is  based  on  a  direct 
analogy  to  this  method,  a  review  is  necessary  of  the  work  done  by  Li 
and  Geiger  in  Ref.  4.  This  method,  as  opposed  to  other  more  exact 
methods,  lends  itself  to  simple  calculation. 

The  basic  equations  of  inviscid  fluid  motion  are  the  starting  point 
of  the  analysis.  Certain  assumptions  are  then  made  wliich  allow  their 
simplification  to  a  degree  where  a  closed  form  solution  is  possible. 

The  solution  permits  the  specification  of  the  flow  parameters  through¬ 
out  the  stagnation  region. 

The  momentun  equations,  written  in  terms  of  orthogonal  curvi¬ 
linear  coordinates,  for  invis.,ir,  iTow  are 

U  ♦  <1  +Ky)v  t  Kuv+  i  =0  (13) 

u  ♦  (1  ,Ky)v  -  Ku‘  .  i  =0  (14) 

The  continuity  equation,  reduced  to  two  dimensional  flow,  is 

^  (pu)  +  ^  [{1  +  Ky)pv]  =  0  (15) 

TR  448  15 


CONFIDENTIAL 


CONFIDENTIAL 


The  curvilinear  coordinate  syatem  ia  shown  below  where  the  origin  i» 
taken  at  the  stagnation  point  withy  measured  normal  to  the  body  and  x 
along^the  body.  The  parameter  K  which  may  be  a  function  of  x,  is  the 
local  curvature 


of  the  body  (K=  — d0/dx=  -  I/Rb^’ 

In  the  stagnation  region  of  the  body  Eqs.  (13),  (14)  and  (15)  may 
be  simplified  by  use  of  the  following  conditions: 

1.  In  the  stagnation  region  u/u^.,  v/u^  =  0  (k)  where  k  (density 
ratio  across  shock)  is  very  small-  Thus  higher  order  terms 
Kuv,  Ku*,  etc-,  can  be  neglected. 

2.  Between  the  shock  and  body  Ky  =  0(5/Rb]  and  can  be 
neglected  for  6/R3  «  1. 

3.  For  M  >3  the  Mach  number  immediately  behind  the  shock 

rapidly  approaches  -/(y  -  Thus,  the  flow  field  in  the 
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stagnation  region  nnay  be  assumed  incompressible  p  “  Pz  =  P^/k- 
These  simplifications  reduce  the  basic  equations,  Eqs.  (13),  (14)  and  (15)  to 

u  (8  u/8  x)  +  V  (8  u/8  y)  =  (  -  k/p*  )  O  p/9  x)  (1^6) 

u  (8  v/8  x)  +  V  (8  v/8  y)  =  {-  k/poo  )  O  p/»  y) 

(8u/8x)  +  (0  v/8y)  =  0  (18) 

A  solution  of  these  differential  equations  must  be  found  subject 

to  the  following  boundary  conditions: 

at  the  body  y  =  0,  v  =  0,  and  at  a  point  A  just  behind  the  shock  wave 
y  =  6  “a  "  ‘^A  ~ 

V.  =  ~qA  si*^  (8  -  (^^) 

A  A 

=  N/(kU(6  sin  /3)*  +  (ut*  cos  W 

P^  =  Poo  +  PooUoo*  (1  -  k)  sin^3 

Also,  at  the  stagnation  point  (x  =  0,  y  =  0),  u  =  v  =  0. 

Near  the  stagnation  point  it  is  assumed  that 

u  =  xf  (y)  =  X  (df/d  y)  (20) 

It  is  seen  that  Eq.(18)  is  satisf.ed  if 

V  =  - f (y)  (^^ ) 

Eliminat-on  of  the  pressure  from  Eqs.  (16)  and  (17)  results  in 
the  equation 

u  Ou/8  x)  +  v(8u/8y)  =  0  (22) 


where  the  vorticity,  «,  is  written 


w  =  (8  v/8  x)  -  (8  u/8  y) 
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Equations  (20),  (21)  and  (22)  yield  a  total  differential  equation 
for  f  (y)  which  is 

f  f  -  £  £'"  =0  123) 

This  is  eUi  ordinary  third  order  differential  equation  which  has  the 
general  solution 

1:  £(y)  =  Ao  exp  (A^  y)  +  A,  exp  (-Aj  y) 

VV 

r  where  the  constants  of  integration  Aq  ,  A|  and  Aj  are  determined  by 

the  boundary  conditions.  An  evaluation  of  the  constants  in  the  region 
near  the  stagnation  point  (i.  e.,  *♦  0  and  A  -►  C  on  sketch)  allows 

limiting  values  of  the  angles  to  be  taken  and  6  -♦  S*. 

Substitution  of  the  constants  into  Eqs.  (20),  (21)  and  (22)  yields 
the  following  results: 


u/uflQ  =  (-x/Rg)  N/k(2  -  k)  cosh  ^(y/S'i*)  sinh  *  [(1  -k)/'dk(2  -  k)]-^ 

(24)  • 

v/vigg  =  -k  [’vyk(2  -  ki  /(I  -  k)]  8inh-^(y/6'«‘)  sinh  ' 

[(1  -  k)/NfkC2  ^k)  ])-  (25) 

a>  =(x/6’i-)  'Vk(2  -  k)  sinh"*  sinh 

-|'(y/6’i‘)  sinh  *  [{1  -  k)/  Vk(2  -  k) 

It  is  also  found  as  part  of  the  solution  that 


b’*  -  [{-kRg)/{l  -  k)]  sinh  *  [(1  -  k)/'s/k(2  -  k)  ] 

and 


(26) 


Cp  =  (2  -  k)  [1  -  (x/Rg)*]  (27)  r, 

i 

Equations  (24),  (25),  (26)  and  (27)  yield  simple  determination  of  any  3 

of  the  flow  variables  in  the  stagnation  region  of  a  blunt  body  in  hyper-  ]  i 

sonic  flov/.  f*  i 


Although  the  assumptions  made  hold  only  near  the  stagnation 
point  in  hypersonic  flow  (i.e.,  K  «  1),  the  theory  gives  a  good 
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approximation  to  pressure  and  velocity  data  for  stream  Mach  numbers 
as  low  as  2.0.  Also,  the  approximate  theoretical  limit  of  the  theory  is 
at  an  angular  position  4)  =  6*;  however,  good  results  ar-e  obtained  for 
values  of  ijj  up  to  about  30*. 

2.  Primary  Flow  Field 

The  flow  field  of  the  primary  stream  may  be  considered 
directly  analogous  to  the  blunt  body  problem  for  injection  normal  to  a 
supersonic  stream.  This  sinalogy  is  based  on  the  assumption  that  no 
mixing  occurs  between  the  primary  and  secondary  streams.  In  this 
case  there  is  some  boundary  separating  the  streartis.  Restricting  the 
problem  to  two  dimensions,  the  boundary  may  be  considered  analogous 
to  the  wall  of  a  cylinder.  The  primary  flow  field  is  then  analogous  to 
the  flow  past  a  cylinder.  A  comparison  of  the  sketches  below  makes 
the  analogy  clear. 


BLUNT  BODY 


SECONDARY  INJECTION 


Boundofy 


7777777 


Secondory  Flow 


jectlon 
to  Woll 


Nor  mol 


Figure  3  is  a  schlieren  photograph  of  injection  into  the  supersonic 
portion  of  a  nozzle  at  an  angle  60"  upstream  of  the  normal  to  the  wall. 
Although  the  photograph  does  not  seem  to  substantiate  the  sketch  shown 
above,  an  explanation  is  apparent.  The  analogy  assumes  inviscid  flow 
whereas  in  reality  the  viscous  fluid  boundary  layer  permits  disturbances 
to  be  transmitted  upstream  along  the  wall.  The  disturbances  separate 
the  boundary  layer  from  the  surface  causing  an  oblique  shock  to  emanate 
from  the  point  where  separation  occurs. 
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Closer  examination  of  the  photograph,  ho'wever,  r^ve-'  a 
curved  shock  formation  (arrow  in  Fig.  3)  which  woulf*  oo  preifnt  lor 
aljluht  body  in  low  Mach  number  supersonic  flew.  Tiierefor.",  v  - 

may  be  regardedas  lending  some  support  to  the  blunt  '  dy  i  naloj  y. 

Equations  (24),.  (25),  (26)  and  (27)  of  th»*  j/revioue  section  yield 
values  for  all  the  primary  flow  variables  in  the  regif>n  corr  spondin,/,  to 
the  stagnation  region  of  the  blunt  body  once  the  free  s^'  eair.  Mac:, 
number  and  shape  of  the  secondary  flow  region  are  specifiei'.  Ti.e 
criteria  for  this  specification  are  discussed  below  under  th  i  headi;  g 
of  Combined  Flow  Field. 

3.  Secondary  Flow  Field  j 

The  secondary  flow  fieV  is  contained  within  the  ref^ion  occ”  li . 
by  the  body  in  the  blunt  body  problem  and  is,  in  fact,  the  l/lunt  bod^'  of 
the  analogy.  Thus  the  problem  consists  of  trying  to  find  a  flow  fi^lc 
occupy  the  region  marked  secondary  flow  in  the  sketch  on  page  l*7j..'  On’.y 
one  possible  flow  field  will  be  discussed  . 

Consider  the  following  sketch  of  a  ’  jconoai  ,  ..ii  "  ■. 

Boundary  Between  F'ri'nory 


The  secondary  flow  is  injected  through  slot  of  wid.n  i  with  velocity  Ujng 

pressure  ,  and  density  p„  •  It  is  a.^su:..  that  the  secondary  flow 
s  s 

Vv 
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A  further  assumption  is  now  made  that  the  secondary  stream 
is  a  perfect  gas  in  isentropic  flow.  This  permits  the  use  of  the  pressure 
density  relation  j  - 

p  =  const  (36) 

With  this  assumption  the  pressure  and  density  may  be  eliminated 
from  Eqs.  (28)  and  (29)  by  cross  differentiation  to  yield: 


whereK  is  considered  to  be  a  function  of  x.  Combining  Eqs.  (33),  (34) 
and  (37)  yields  a  total  differential  equation  for  K  and  f(x)  which  is 

y(ff'"  -  f  f)  +K'  f2  +Kff'  =  0  (38) 

where  the  primes  denote  differentiation  with  respect  to  x.  It  should 
be  noted  that  Eq.  (38)  is  similar  to  Eq.  (23)  obtained  by  Li  and  Gicger 
except  for  the  terms  containing  K  and  its  derivative.  These  terms 
arise  from  nhe  fact  that  v/u^Og  and  u/u^Og  are  not  necess.irily  small 
compared  to  unity.  In  order  to  solve  Eq.  (38)  various  K  functions  have 
to  be  aswurned  and  the  corresponding  f  functions  determined  numerically. 

4.  Combined  Flow  Field 

The  remaining  bound.iry  conditions  on  the  secondary  flow 
field  are  imposed  when  the  flows  are  combined.  Since  the  steady 
state  problem  must  be  in  equilibrium,  it  is  necessary  that  the  static 
pressures  of  the  primary  and  secondary  streams  match  across  the 
boundary.  For  any  particular  K (x)  the  pressure  on  the  boundary 
exerted  by  the  primary  flow  is  specified  when  the  values  cf  pao,  pooi 
and  UjiQ  are  known  in  the  primary  flow.  For  the  particular  K  (x) 
chosen  the  function  f  (x)  is  determined  and  therefore  the  pressure  of 
the  secondary  flow  on  the  boundary. 
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fills  only  n.  .  the  region  between  the  wall  and  the  boundary 

o  ■’A  -'y  and  secondary  flows-  The  remainder  of  the 

region  is^a  zone.  It  i«  also  assumed  -  — 

ky  =  0(./H,;'  r:  •  .  be  leglected.  However,  u/u^q  and  v/v^q  may  hot 
be  neglecieoi  'f.  <  e  s.  ^  ondary  flow  as  they  were  in  the  primary  flow. 

Ht)'.;  ;13  ;  (14)  and  (15)  are  simplified  by  the  assumptions 

made  abo'  ;aay  .  i  written  as  follows: 


1  .''v./O'  +  V  (du/8y)  +  Ku  V  =  - (28) 

+  v(8v/9y)  -Ku*  =  -  ^  (29) 

-  (pu)  0  +  9  (pv)/0  y  =  0  (30) 


Thi  o  dUr<;.i  •‘■jal  oquat  ons  must  be  solved  to  yield  solutions  subject 
to  cert  'i  -  b  nmdax'y  con-  itions.  Two  of  the  required  boundary  conditions 
arise  from  ti^  geometry  of  the  sketch  and  are  as  follows. 

^^\e  tai^ency  condition  at  the  boundary  between  the  primary 
a.  •’  8ec|5ndury»strei.m8 

k=  0  (31) 


I  ti' 

and  the  cond/tions  at  the  injection  port 


x=0,  \i  =  Uoo  d<y<0 


(32) 


The  I  (■luaiiiing  botmdary  conditions  are  discussed  when  the  flow  fields 
are  corriutued  in  the  next  section. 

la  the  thin  region  of  secondary  flow,  it  may  be  assumed  that 

u  =  fix)  (33) 


V  =  -  yf  (x)  =  -  y(df/dx)  (34) 

.'he  continuity  equation,  Eq .  (30)  is  satisfied  if 

p  =  Poo  +  yf 

I'he  choice  of  thx  particular  form  of  v  satisfies  Eq.  (31). 
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A  further  assumption  is  now  made  that  the  secondary  stream 
is  a  perfect  gas  in  isentropic  flow.  This  permits  the  use  of  the  pressure- 
ienrity  relatioh  - 

p  =  const  (36) 

With  this  assumption  the  pressure  and  density  may  be  eliminated 
from  Eqs.  (28)  and  (29)  by  cross  differentiation  to  yield: 


whereK  is  considered  to  be  a  function  of  x.  Combining  Eqs.  (33),  (34) 
and  (37)  yields  a  total  differential  equation  for  K  and  f  (x)  which  is 

y(ff"i  _  f.  f.i)  +K'  f2  +Kff'  =  0  (38) 

where  the  primes  denote  differentiation  with  respect  to  x.  It  should 
be  noted  that  Eq.  (38)  is  similar  to  Eq.  (23)  obtained  by  Li  and  Gieger 
except  for  the  terms  containing  Kand  its  derivative.  These  terms 
arise  from  the  fact  that  and  u/u^,^  are  not  necess.arily  small 

compared  t<;  unity.  In  order  to  solve  Eq.  (38)  various  K  functions  have 
to  be  asHurned  and  the  corresponding  f  functions  determined  numerically. 

4.  Combined  Flow  Field 

The  remaining  bound-iry  conditions  on  the  secondary  flow 
field  are  imposed  when  the  flows  are  combined.  Since  the  steady 
state  problem  must  be  in  equilibrium,  it  is  necessary  that  the  static 
pressures  of  the  primary  and  secondary  streams  match  across  the 
boundary.  For  any  particular  K (x)  the  pressure  on  the  boundary 
exerted  by  the  primary  flow  is  specified  when  the  values  cf  p*,  pgo, 
and  are  known  in  the  primary  flow.  For  the  particular  K  (x) 
chosen  the  function  f(x)  is  determined  and  therefore  the  pressure  of 
the  secondary  flow  on  the  boundary. 
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Since  a  nvimerical  solution  was  not  obtained,  it  is  not  known 
whether  an  exact  pressure  match  across  the  boundary  is  possible  by 
varying  K(x).  At  best  it  would  seem  that  an  approximate  pressure 
match  over  some  portion  of  the  boundary  would  have  to  be  met  by 
this  approach. 

The  secondary  stream  pressure  on  the  boundary  as  a  function 
of  X  is  determined  by  the  relation 


Ps  =  R»8 


(39) 


where 


r  -  V  -  R  -  V 

r-y  K 


thus 


Ps 


Pcos  + 


.1: 


pg  (Ky  +  1) 
K 


u*  d  y 


(40) 


The  lower  limit  of  the  integral  in  Eq.  (40)  is  the  value  of  y  on  the 
interface  between  the  secondary  flow  and  the  separated  region  for  a 
particular  value  of  x.  The  interface  is  the  path  followed  by  a  particle 
being  injected  from  the  edge  of  the  slot  y  =  -  d. 

The  analogy  is  limited  to  small  values  of  the  angle  due  to  the 
fact  that  the  analysis  of  Li  and  Geiger  is  valid  only  near  the  stagjiatlon 
point  of  a  blunt  body.  However,  the  analogy  might  be  useful  as  a 
means  of  examining  a  portion  of  the  flow  field  generated  by  secondary 
injection  normal  to  a  supersonic  stream. 


D.  FURTHER  ANAL.YTICAL  RESULTS-NSL  FLOV,'  MODELS 

The  Phase  I  data  corroborated  the  sheet  flow  and  secondary 
expansion  models  to  an  encouraging  extent.  However,  all  of  the  data 
corresponds  to  the  case  of  .secondary  air  injection  into  a  primary  strean 
wliich  is  also  air.  In  addition,  the  ratio  of  primary  to  secondary  stag¬ 
nation  temperature  was  approximately  unity.  The  range  of  the  experi¬ 
mental  data  obviously  does  not  provide  exhaustive  coverage  of  the 
\-ariables  but  it  should  also  be  noted  tnat  this  was  not  the  intention  of 
the  investigation.  Nevertheless,  in  spite  of  a  lack  of  data,  some 
attention  was  given  to  the  effects  of  differences  in  gas  and  temperature 
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between  the  primary  and  secondary  strearr. 8.  This  work  is  based  on 
calculations  in  which  temperatui'e  and  gas  effects  were  studied  to 
— determine^heir  effect  on  flow  throttling. 

1.  Temperature  Effect 

The  effect  of  temperature  has  been  studied  through  the  use  of 
the  ratio  of  secondary-to-primary  stagnation  temperatures.  It  has 
been  found  that  this  temperature  ratio  plays  an  important  part  in  the 
effectiveness  of  seconv..ary  injection  as  a  means  of  flow  control.  The 
following  discussion  concerns  the  stagnation  temperature  ratio  effect 
on  flow  throttling  performance  as  predicted  by  the  sheet  flow  and 
secondary  expansion  models.  It  should  be  noted  that  nn  change  has 
been  made  in  the  analyses  to  account  for  heat  transfer  between  the 
primary  and  secondary  streams.  In  view  of  the  nature  of  the  flow 
models,  in  which  the  two  streams  flow  without  mixing,  this  seems  a 
satisfactory  approach  for  a  first  approximation. 

The  continuity  and  momentum  equations  for  the  sheet  flow  and 
secondary  expansion  models  are  given  by  Eqs.  (8),  (9),  (10)  and  (12). 

The  stagnation  temperature  ratio  does  not  appear  in  any  of  these  equation 
Thus,  flow  throttling  performance  is  specified  for  a  given  secondary- 
to-primary  stagnation  pressure  ratio.  This  means  that  the  temperature 
ratio  effect  is  felt  only  through  the  expressinn  for  the  secondary  flow 
ratio.  In  particular,  examination  of  the  secondary  flow  ratio  (Eq.  (11)) 
shows  that  it  varies  inversely  with  the  square  root  of  the  secondary- 
to-primary  stagnation  temperature  ratio.  This  variation  is  illustrated 
by  the  results  shown  in  Fig.  4  where  flow  throttling  performance  is 
given  for  a  typical  configuration.  Although  this  is  the  somewhat  ideal 
case  of  air  injection  into  an  air  stream,  the  trend  is  nevertheless 
clear  regarding  the  decreasing  effectiveness  of  injection  as  the  stag¬ 
nation  temperature  ratio  decreases.  Furthermore,  at  a  temperature 
ratio  of  0.  2  the  level  of  constant  total  flow  is  being  approached  wherein 
the  thrust  remains  approximately  constant.  Morover,  this  flow 
throttling  behavior  is  discouraging  in  view  of  the  fact  that  in  a  real 
application  a  practical  value  of  the  temperature  ratio  would  probably 
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be  less  than  unity. 

At  present  there  is  no  flow  throttling  data  available  for  the  case 
of  injection  at  the  throat  in  which  the  stagnation  temperature  ratio  differs 
significantly  fiom  unity.  Rocketdyne  reported  in  Ref.  5^at  secondary 
injection  just  upstream  of  the  throat  at  a  stagnation  temperature  ratio  of 
0. 15  did  not  throttle  the  primary  flow.  This  seems  to  disagree  somewhat 
(although  not  greatly)  with  the  trend  predicted  by  the  NSL  flow  models. 
However,  the  low  secondary  flow  rates  (wg/wp^  up  to  a  maximum  of 
0. 04)  and  subsonic  injection  associated  with  the  Rocketdyne  tests  are  of 
sufficient  importance  to  account  for  the  discrepancy. 

It  should  be  mentioned  that  UAC  (Ref.  6)  found  the  effectiveness 
of  vector  control  by  injection  to  vary  with  the  square  root  of  the  secondary- 
to-primary  stagnation  temperature  ratio.  Furthermore,  the  form  of  this 
dependence  is  the  same  as  that  of  the  NSL  flow  models.  That  is,  perfor¬ 
mance  (some  measure  of  side  force)  as  a  function  of  secondary  flow 
varies  in  the  same  manner  as  flow  throttling.  Thus,  there  is  some 
support  (admittedly  meager)  for  the  effect  of  temperature  given  by  the 
NSL  models.  However,  until  conclusive  flow  throttling  data  is  available 
the  predictions  of  the  effect  of  temperature  must  be  regarded  as 
unsubstantiated . 

2.  Injection  Gas  Effect 

The  generalized  form  of  the  continuity,  momentum,  and  secondary 
flow  equations  for  the  sheet  flow  and  secondary  expansion  models  permits  the 
variation  of  the  gases  in  the  primary  and  secondary  streams.  Examina¬ 
tion  of  these  equations  reveals  chat  a  low  ratio  of  specific  heats  and  low 
molecular  weight  are  desirable  properties  for  the  secondary  gas.  That 
this  is  indeed  the  case  is  shown  by  the  flow  throttling  results  in  Fig.  5. 
Helium  injection  into  air  is  shown  for  a  typical  configuration’^.  Several 
stagnation  temperature  ratios  are  considered  and  the  results  for  air 
injection  are  shown  for  compari  son.  Helium  injection’^’!*  obviously  makes 

‘^The  curves  shown  for  helium,  like  those  for  air,  start  at  the  point  where 
seccnc  ry  injection  becomes  sonic,  as  explained  in  the  last  paragraph 
r'-  pa-e  8. 

'^’■‘Although  helium  has  a  high  y,  its  low  molecular  weight  makes  it  a 
good  injectant  since  the  effect  of  molecular  weight  is  stronger  than  that 
of  y . 
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i\  difference  in  flow  throttling  performance  being  a  good  deal  better 
than  air  injection.  In  spite  of  such  encouraging  results,  there  is,  once 
again,  no  data  to  support  or  deny  these  predictions. 

3.  Combustion  Effect 

The  effect  of  combustion  has  been  very  crudely  estimated  for 
the  sheet  flow  and  secondary  expansion  models.  This  was  done  by 
substituting  real  engine  values  of  the  specific  heat  ratio  and  molecular 
weight  (gas  constant)  for  the  primary  stream  in  the  appropriate  equations. 
The  results  are  shown  in  Fig.  6  where  flow  throttling  is  given  for  several 
cases’*'.  Helium  injection  into  hydrogen  peroxide  and  hydrazine  or  liquid 
oxygen  and  gasoline  is  not  as  good  as  helium  injection  into  air.  A 
stagnation  temperature  ratio  of  0.2  was  chosen  since  this  corresponds 
to  a  case  where  helium  at  1200 *R  might  be  injected  into  a  6000'R  chamber. 
In  spite  of  the  decrease  in  effectiveness  due  to  the  combustion  gas 
parameters  (and  low  injection  temperature),  helium  injection  still  shows 
some  potential  for  flow  throttling.  A  mentioned  before,  no  data  is 
available  for  comparison. 

E,  PHASE  II  --  TWO-DIMENSIONAL  NOZZLE  TESTS 

The  Phase  n  portion  of  the  NSL  program  was  a  study  of  secondary 
injection  into  two-dimensional  nozzles.  One  objective  of  the  tesla  in¬ 
volved  was  to  determine  the  effect  on  flow  throttling  of  various  param¬ 
eters  whose  importance  had  not  been  thoroughly  established  by  the 
Phase  I  experiments.  These  parameters  were  injection  angle,  nozzle 
throat  region  pressure  gradient,  and  nuzzle  pressure  ratio.  The  other 
objective  was  to  gain  some  understanding  of  the  basic  nature  of  the  flow 
throttling  process  through  schlicren  photography.  Schlieren  flow  obser¬ 
vations  were  made  through  glass  windows  wliich  formed  the  major 
portion  of  the  nozzle  sidewalls. 

See  first  footnote  on  page  25. 
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The  first  part  of  this  section  describes  the  model  and  test  con¬ 
ditions.  The  remainder  is  a  presentation  and  discussion  of  the  data. 

The  vector  control  aspect  of  the  Phase  II  study  is  presented  in  a 
subsequent  section. 

1.  Model  Details  and  Test  Conditions 

a.  Air  Supply  System 

The  schematic  diagram  of  Fig.  7  shows  the  supply  system  for 
both  the  primary  and  secondary  air.  Chicago  Pneumatic  Compressors 
pumped  the  imary  air  through  a  drier  into  three  storage  tanks  denoted 
in  the  diagram  as  receivers.  A  main  control  valve  throttled  the  air 
as  it  entered  a  surge  tank  before  passing  to  the  primary  flowmeters. 

The  flowmeters  were  used  to  measure  the  primary  mass  flow  of  air 
through  the  model.  From  the  flowmeters  the  air  flowed  to  a  settling 
chamber  at  the  model.  Auxiliary  vacuum  pumps  exhausted  the  nozzle 
flow  after  it  passed  through  a  diffuser.  A  thermocouple  and  pressure 
gage  located  in  the  settling  chamber  were  used  to  measure  the  primary 
stream  stagnation  conditions.  The  main  control  valve  was  operated  so 
as  to  maintain  a  constant  -itagnation  pressure  irrespective  of  the  amount 
of  primary  flow  throttling.  Stagnation  pressures  of  40,  50,  60  and  70  psia 
were  used  during  these  tests. 

Secondary  air  for  these  experiments  was  obtained  from  an 
Inger soli- Rand  compressor  which  pumped  air  into  SOOOpsia  storage  tanks. 
Tliis  system  supplies  air  at  a  working  pressure  of  600  psia.  The  mass 
flow  rate  of  secondary  air  was  measured  by  an  orifice  plate  in  combin¬ 
ation  with  pressure  measurements  upstream  and  downstream  of  the 
orifice.  A  manually  oper.ated  control  valve  was  used  to  regulate  the 
secondary  air  flow  to  obtain  a  desired  stagnation  pressure.  A  3000  psia 
Heise  pressure  gage  was  used  to  measure  the  secondary  stagnation 
pressure  in  the  model  settling  chambers.  The  secondary  stagnation 
temperature  was  obtained  by^  thermocouple  as  in  the  primary  flow. 

b.  Model  Arrangement 

'I  he  model  was  arranged  approximately  as  it  is  shown  in  the 
schematic  diagram  of  Fig.  7.  Since  no  forces  were  being  measured 
directly,  the  model  was  mounted  rigidly  in  place. 
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c.  Nozzle  Design 

Three  nozzles  were  tested  in  the  Phase  II  experiments.  The 
first^ne  tested^was-conical,  as  siiowninPig.  8,  and.it  was  loundlhat 
a  shock  pattern  originated  f rom  overexpansion  at  the  throal^  Since  this 
shock  pattern  interferred  with  flow  observations,  the  remaining  two 
nozzles  were  designed  for  essentially  shock-free  flow.  These  nozzles 
were  designed  mainly  to  have  widely  different  throat-region  pressure 
gradients.  The  "high"  and  "low"  throat-region  pressure  gradient 
nozzle  contours  are  shown  in  Figs.  9  and  10,  respectively.  The  high 
pressure  gradient  is  four  times  greater  than  the  low  one  and,  inciden¬ 
tally,  the  pressure  gradient  of  the  nozzle  in  Fig.  Sis  5.5  times  that  of 
Fig.  10  Uow  one).  The  nozzle  in  Fig.  10  is  a  portion  of  the  contour  of 
the' Mach  number  3. 5  main  NSL.  wind  tunnel  nozzle  blocks  (Ref.  7) 
whereas  the  nozzle  of  Fig.  9  was  designed  according  to  the  method  of 
Ref.  8.  Aside  from  the  throat-region  pressure  gradient,  the  only 
difference  in  nozzles  is  the  exit-to-throat  area  ratio  (which  determines 
the  exit  Mach  number).  The  exit-to-throat  area  ratios  (and  exit  Mach 
numbers)  for  the  nozzles  in  Figs.  8,  9,  and  10  are  5.38(M=  3.25), 

4.  34  (M  =  3.  02),  and  2. 89  (M  =  2.  60),  respectively.  All  three  nozzles 
were  designed  for  3.  mass  flow  of  2  pounds  per  second  at  a  stagnation 
pressure  of  70  psia  and  a  stagnation  temperature  of  530' "R. 

The  contoured  nozzles  were  designed  toaccommodate  interchangeable 
injection-section  inse-vts  at  the  throat  and  in  the  supersonic  portion  of 
the  nozzle  in  Fig.  10.  Injection  angle  was  varied  through  the  use  of 
different  sets  of  insert;;  all  of  the  important  features  of  which  are  shown 
by  the  one  in  Fig.  1 1 . 

Each  nozzle  was  instrumented  with  a  number  of  static  pressure 
taps  along  the  middle  of  the  nozzle  surface.  The  pressui’e  tap  locations 
are  shown  approximately  by  the  ticks  on  the  nozzle  contours  in  Figs.  8,  9, 
and  1C.  Other  pressure  taps  were  used  to  measure  the  stagnation 
pressure  upstream  and  diffuser  pressure  downstream  of  the  nozzle. 

The  windows,  which  were  an  integral  part  of  the  nozzle  side- 
walls,  were  1.  25  inches  thick  w'ith  optically  flat  ground  surfaces.  The 
window  diameter  of  8  inches  (see  Figs.  8,  9,  and  10)  permitted  flow 
observation  through  most  of  the  nozzle. 
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d.  Injection  Configurations 

The  nozzle  and  injection  configurations  tested  are  listed  in  the 
table  below.  In  all  cases  the  nozzle  throat  was  1. 5"  in  height 
in  width.  The  injection  slot  extended  the  width  of  the  nozzle  and  was 
0. 08"  in  the  nozzle-axis  direction  measured  normal  to  the  injection 
slot  axis  (see  Fig.  11). 


Configuration’^' 


Nozzle 


N 

45  °  convergent 

25 ”  divergent 

60 

40, 

50. 

60, 

70 

P 

10”  convergent 

25  °  divergent 

60 

70 

Q 

high  pressure 
gradient  contoured 

60 

40, 

50. 

60, 

70 

R 

same  as  Q 

30 

40. 

50. 

60, 

70 

s 

low  pressure 
gradient  contoured 

60 

50, 

70 

T 

same  as  S 

45 

50, 

70 

U 

same  as  S 

30 

50, 

70 

V 

same  as  S 

75 

50. 

60 

Injection  at  the  inflection  point  in  the  divergent 
portion  of  the  no-.’zle 


W  '’ame  as  S  45  40,  50,  60 

Injection  at  the  inflection  point  in  the  divergent 
portion  of  the  nozzle 


’■■■  Unless  otherwise  specified  all  injection  locationE  -re  at  the  nozzle 
throat  The  injection  angle  is  the  angle  included  between  the  primary 
nozzle  axis  and  the  axis  of  injection.  The  configurations  listed  are 
ail  two  dimensional. 

Asymmetric  injection  data  was  obtained  for  each  of  the  configurations 
listed  by  simply  cvitting  off  t.hc  air  supply  to  one  of  the  injection  slots.  The 
subscript  s  is  used  to  denote  an  asym.metric  configuration  (for  example  Vg). 
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2.  Influence  of  Test  Parameters 

The  primary  objective  of  the  two-dimensional  tests  was,  as 
mentioned  before,  to  determine  the  effecla  of  certain  parameters  on 
flow  throttling  performance.  These  parameters  were  injection  angle, 
throat-region  pressure  gradient,  and  injection  station  and  their  influence 
on  flow  throttling  is  discussed  in  this  section. 

Although  the  Phase  I  test  results  showed  that  nozzle  pressure 
ratio  was  not  an  important  parameter,  a  check  on  this  was  made  in  the 
Phase  II  tests.  As  shown  in  the  configuration  list  of  the  previous  section 
several  primary  stagnation  pressure  levels  were  tested  for  most  models. 
The  data  shown  in  Figs.  12,  13  and  14  is  typical  of  all  of  the  data  as 
regards  the  effect  of  primary  stagnation  pressure.  The  flow  throttling 
results  for  configurations  Q,  R.and  S  are  plotted  in  the  form  of  primary 
and  secondary  flow  ratios  normalized  by  the  appropriate  value  of 
unrestricted  primary  flow.  It  is  evident  from  the  data  in  Figs.  12 
and  13  that  there  is  no  effect  of  prim^lry  stagnation  pressure  on  the 
flow  throttling  of  configurations  Q  aiKl  R.  There  is  a  slight  effect 
observable  in  Fig.  14  where  the  data  of  configuration  S  indicates  that 
flow  throttling  is  increased  by  a  decrease  in  po  from  70  to  50  psia. 

The  effect  is  small,  however,  and  agrees  with  the  trend  reported  in 
Ref.  9  which  is  discussed  further  in  Ref.  1.  It  is  sufficient  to  state 
here  that  the  primary  stagnation  pressure  does  not  affect  flow  throttling 
if  the  overall  nozzle  pressure-ratio  and/or  the  throat-region  pressure 
gradient  are  .sufficiently  large. 

a.  Injection  Angle 

The  effect  of  injection  angle  (^)  on  flow  throttling  was  measured 
with  the  high  and  low  pressure  gradient  nozzles.  The  flow  throttling 
data  for  corresponding  configurations  is  typified  by  that  shown  in 
Figs.  15  and  16.  The  data  in  these  figures  corresponds  to  injection  at 
the  nozzle  throat.  Flow  t.nrottling  increases  with  decreasing  injection 
angle  and  is  evidently  independent  of  nozzle  throat-region  pressure 

The  experimental  uncertainties  in  the  flow  measurements  of  the 
Phase  II  tests  are  ±0.  01  in  Wp/wp^^  and  ±0.  005  in  Wg/wp^  . 
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gradient.  In  addition,  the  data  in  Fig.  16  indicates  that  the  increase  in 
throttling  becomes  smaller  as  injection  angle  decreases.  Both  of  these 

trends  substantiate  qualitatively  the  analytical  results. _ In^lhe^analytical 

models  decreasing  injection  angle  means  that  the  axial  corhponeht  of 
injection  momentum  is  increasing  as  indicated  by  the  cos  ^  in  Eqs.  6, 

10. and  12.  Increases  in  this  injection  stream  momentum  decrease 
the  primary  flow  ratio  (increase  flow  throttling).  Furthermore,  since 
this  momentum  varies  as  cos  (j>  flow  throttling  increases  should  become 
smaller  as  iji  becomes  smaller.  Quantitative  support  of  the  theoretical 
^  effect  is  not  possible  in  view  of  the  uncertainty  in  the  data;  however, 
it  is  believed  that  the  data  does  indicate  a  physically  correct  analytical 
representation. 

Injection  into  a  nozzle  at  a  station  where  the  flow  is  supersonic 
will  obviously  not  throttle  the  flow  until  the  secondary  flow  is  large.  This 
is  60  not  only  because  the  primary  stream  has  more  momentum  at  a 
supersonic  station  but  also  because  the  secondary  stream  expands  down¬ 
stream  more  rapidly  due  to  the  increasing  area  in  that  direction.  The 
data  shown  in  Fig-  17for  supersonic  stationinjectioninto  the  low  pressure 
gradient  nozzle  shows  the  large  secondary  flow  required  for  choking.  It 
is  perhaps  obvious  that  this  data  also  points  out  that  flow  throttling  should 
be  controlled  at  the  nozzle  throat.  Also' apparent  in  Fig.  17  is  the  fact 
that  injection  angle  partially  determines  the  secondary  flow  rate  for 
choking.  Finally,  once  the  nozzle  is  choked  primary  flow  throttling 
varies  with  secondary  flow  as  though  injection  were  at.  the  throat  (which, 
in  effect,  it  is). 

b.  Pressure  Gradient 


Nozzle  throat- region  pressure  gradient  was  investigated  as  a 
factor  which  might  partially  determine  the  secondary,  flow  necessary  to 
separate  the  primary  flow  and  achieve  a  given  amount  of  throttling. 
Typical  pressure  distributions  for  the  high  and  low  pressure  gradient 
nozzles  are  shown  in  Figs.  18  and  19.  As  mentioned  before,  the  throat- 
region  pressure  gradients  are  in  the  ratio  of  about  4  to  1 .  The  pressure 
distributions  which  correspond  to  flow  wdth  secondary  injection  are 
discussed  in  connection  witli  flow  observations  in  a  later  section- 
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The  effect  of  pressure  gradient  on  flow  throttling  is  illustrated 
in  Figs.  20  and  21  for  the  injection  angles  of  30°  and  60°,  The  effect 'is 
appareixtly  small nnd  independent  of  injection  angle.  A  decrease-in 
threat-region  pressure  gradient  produces  an  increase  in  flow  throttling 
and  this  change  in  throttling  is  very  slightly  attenuated  by  increasing  , 
stagnation  pressure  (po). 

The  small  increase  in  flow  throttling  which  results  from  a  reduc¬ 
tion  in  pressure  gradient  may  be  explained  by  the  fact  that  the  primary 
'  stream  is  actually  not  sensitive  to  nozzle  contour  once  there  is  a  large 
secondary  flow.  That  is,  large  secondary  flow  generates  an  effective 
throat  which  throttles  the  primary  flow  and  changes  its  path  considerably. 
Since  the  secondary  flow  thus  provides  an  effective  or  "fluid"  nozzle, 
the  primary  flow  does  not  feel  the  wall  directly.  Therefore,  the 
"pressure  gradient"  influence  is  actually  a  geometry  one  in  which  the 
effective  ("fluid")  nozzle  is  partially  determined  by  the  physical  nozzle 


throat. 


c.  Injection  Stat;  on 


The  influence  of  injection  station  was  mentioned  in  part  above 
but  will  be  reviewed  here  for  chirity  and  completeness.  Tlie  data  in 
Fig.  17  shows  that  primary  flow  throttling  does  not  take  place  until 
the  nozzle  is  choked  by  secondary  injection.  This  indicates  further 
that  the  nozzle  throat  i.s  the  station  at  which  throttling  should  be  most 
efficient  (require  the  least  secondary  flow).  The  data  also  indicates 
that  tlirottling  is  the  same  function  of  secondary  flow  once  the  nozzle 
is  choked.  This  again  indicates  injection  should  be  at  the  throat 
since  that  avoids  the  need  for.  the  additional  secondary  flow  required 
for  choking. 

d.  Thrust  Control 

In  a  conventional  nozzle  the  thrust  is  directly  proportional  to 
the  product  of  mass  flow  and  exit  velocity.  The  thrust  of  a  nozzle  being 
throttled  by  secondary  injection  has  similarly  been  found  to  be  propor¬ 
tional  to  the  total  (secondary  plus  pr'inary)  flow  (Refs.  1  and  9).  The 
Phase  II  test  data  also  agrees  with  this  result  as  described  on  the 
following  page. 
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The  value  of  a  thrust  parameter,  defined  in  Ref.  1  as 

X 

Fq 

w  +  w 

P  a 

Wr. 

Po 

was  calculated  for  a  large  sample  of  the  Phase  II  tests.  The  thrusts  F 
and  Fo  corresp  nd  to  the  flow  rates  Wp+Wg  and  Wp^ ,  respectively,  and 
were  calculated  from  pressure  integrations  together  with  estimates  of 
axial- component  jet  reaction  forces  (of  the  secondary  flow).  All  of  the 
thrust  parameter  values  fell  between  0.98  and  1.07  and  the  sampling 
included  asymmetric  as  well  as  symmetric  injection.  This  range  agrees 
closely  with  the  corresponding  one  for  the  Phase  I  tests  (Ref.  1)  in 
which  the  values  fell  between  0.96  and  1. 04. 

The  variation  of  the  thrust  ratio*^'  (F/Fo)  and  total  flow  ratio 
(wp+  Wg/wp^)  are  shown  as  a  function  of  secondary  flow  in  Figs.  22  and 
23.  The  data  in  Fig.  22  shows  that  results  are  independent  of  stagnation 
pressure  (po).  an  observation  common  to  all  configurations.  Also 
evident  is  the  fact  that  thrust  is  related  to  total  nozssle  flow  by  a  constant. 
This  i.s  not  so,  however,  for  the  configurations  of  Fig.  23  where  the 
proportionality  varies  with  secondary  flow.  These  results  are  typical 
of  the  other  low  pressure  gradient  nozzle  and  side  injection  configurations 
whereas  those  of  Fig.  22  are  typical  of  the  high  pressure  gradient  nozzle 
configurations.  It  is  apparent  in  Fig.  23  th.it  the  total  flow  decrea.ses 
faster  than  the  thrust  as  secondary  flow  increases,  a  result  which  is 
opposite  to  the  trend  observed  in  the  Phase  I  tests  (Ref.  1).  There  are 
two  obvious  possible  factors  involved  in  this  discrepancy.  First,  the 
Phase  I  tests  concerned  an  axially  symmetric  nozzle  while  the  present 
nozzles  are  two-dimensional.  Second,  the  proportionality  (between 
thrust  and  total  flow)  varies  with  secondary  flow  only  in  configurations 
where  the  nozzle  does  not  expand  to  a  high  exit-to-throat  area  ratio 
except  for  two-dimensional  side  injection.  Tlius,  nozzle  geometry  and 

The  experimental  uncert.iinty  in  F/Fq  is  iO.015 
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expansion  ratio  are  possibly  of  some  importance’!'.  It  is  believed  that  nozzle 
geometry  is  the  key  factor,  at  least  qualitatively,  but  the  quantitative 
nature  of  the^ data  makes  it  inconclusive.  . 

3.  Flow  Observations 

One  of  the  main  purposes  of  using  two'dimensional  nozzles  in  the 
Phase  II  tests  was  to  pernrut  the  use  of  schlieren  photography.  It  was 
felt  that  extensive  flow  observations  would  reveal  the  nature  of  the  throttling 
process.  Although  the  analytical  model  assumptions  were  formulated 
primarily  on  the  basis  of  intuition,  Phase  I  data  indicated  that  the  assumptions 
were  sound-  The  Schlieren  photographs  discussed  in  this  section  further 
afford  visual  substantiation  to  the  throttling  analysis. 

The  flow  patterns  associated  with  various  rates  of  secondary 
injection  arc  shown  in  Fig.  lA.  This  series  of  photographs  concerns  the 
high  pressure  gradient  nozzle  (configuration  Q).  The  wide  disturbance  at 
the  throat  in  Fig.  24a  is  due  to  the  injection  slot  wliile  the  narrow  one 
just  upstream  emanute.s  from  the  insert  joint  (see  Fig.  9).  The  weak 
shock  wave  downstream  of  the  throat  is  caused  by  some  slight  overex- 
pansion  iii  the  throat  region.  In  Fig.  24b  the  over  expansion  must  have 
increased  with  injection  since  the  resultant  shock  wave  i?  stronger.  The 
secondary  flow  evidently  diffuses  without  separation  since  the  overex¬ 
pansion  shock  proceeds  to  the  nozzle  wall.  The  density  gradient  in  the 
secondary  injection  vicinity  is  so  strong  that  flow  details  are  obsctired. 

In  Fig.  24c,  however,  new  flew  features  become  apparent  with  a  higher 
secondary  flow  rate.  The  primary  flow  moves  at  low  velocity  upstream 
01  injection  as  indicated  by  the  faint  circular  waves  which  propagate  from 
the  injection  slots.  A  thick  boundary  layer  is  apparent  just  downstream 
of  injection  and  the  overexpansion  shock  waves  have  increased  in 
intensity.  Figs.  24d  and  24e  make  the  subsonic  region  waves  more 
evident  and  show  even  stronger  overexpansion  shock  waves.  More 
important,  however,  these  photographs  show  clearly  the  features  of  a 
flow  which  is  a  combination  of  mixing,  sheet  flow,  and  secondary  expan- 

It  should  oc  noted  that  in  all  cases  injection  effects  extend  only  a  short 
distance  downstream  of  the  i.ajection  station  as  shown  in  Figs.  18  and  19. 
Thrust  variation  with  injection  in  the  supersonic  portion  of  the  nozzle 
is  thus  tied  strongly  to  geometry  and  expansion  ratio. 
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sion.  The  flow  configuration  appears  to  be  arranged  as  shown  in  the 
sketch  below. 


y- Mixing  RiQtpn 
y/^  y^Sheet  Flow 
y'  y  y-Exponslon  Fon 
y  y  /  y^Region  of  Secondory  Exponslon 

Primary  Flow  - 


Secondary 

Flow 


-Shock  Wove 
from  Overexpansion 


All  of  these  details  are  quite  apparent  in  Figs.  ?Ad  and  24e  and  the  flow 
configuration  as  a  whole  substantiates  the  basic  nature  of  the  throttling 
,1  na  ly  si  s  , 

The  achlieren  photograph  in  Fig.  ZS  shows  a  typical  flow  field 
produced  by  injection  at  the  throat  in  the  low  pressure  gradient  iio/,/,lo. 

The  higher  pressures  of  this  flow  can  be  seen  lu  c  hange  the  angles  of 
the;  overcixpansion  .shock  wave  comp.ired  to  those  in  Fig.  24d.  Higher 
pressure  is  also  responsible  for  the  fact  that  the  expansion  fan  does 
nut  extend  very  far  downstream.  It  is  quite  evident,  however,  that  the 
flow  pattern  is  btisically  *he  same  as  that  in  Fig.  Z4  and  the  configuration 
sketched  above . 

The  scries  of  schliertu  photographs  .shown  in  Fig.  Zb  -corrcspi^nd  to 
iuJei:tion  into  the  supersonic  portion  of  the  low  pressure  gradient  nozile. 
Fluw  scp.iration  upstream  of  the  injection  .-.Idtiuii  is  quite  apparent 
, I  long  with  the  she..;  i.rycr  which  bounds  the  separated  region  .After 
the  pog/pc  ratio  exceeds  unity  the  flow  just  downstream  of  the  injection 
station  becomes  sonic.  Finally,  secondary  flow  chokes  the  nozzle  and 
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moves  the  effective  throat  to  the  injection  location  where  the  flow  is 
being  throttled.  This  final  flow  pattern  (Fig.  26e)  is  easily  recognized 
as^he  same  as  that  in  Figs.  24e  and  25.  - 

Typical  pressure  distributions  rtlong  the  nozzle  contour  are 
shown  in  Figs.  18  and  19  for  the  high  and  low  pressure  gradient  nozzles 
with'secondaTy-injection.  These  pressure  distributions  may  be 
compared  to  the  ochlieren  photographs  in  Figs.  24d  and  25.  This 
comparison  shows  agreement  in  revealing  the  essentially  stagnant 
region  upstream  of  injection  and  low  pressure  expansion  just  downstream. 
Furthermore,  it  is  evident  that  downstream  of  the  over  expansion  shock 
wave  the  pressure  distribution  is  not  very  different  from  that  in  the 
nozzle  without  injection. 

It  should  be  noted  that  an  important  assumption  in  the  theoretical 
models  is  borne  out  by  the  schlieren  photographs  of  Figs.  24,  25,  and  26. 
This  assumption  is  that  the  flow  throttling  process  occurs  in  a  short 
axial  distance  and  thus  may  be  treated  as  taking  place  in  a  constant  area 
duct.  The  flow  observations  as  well  as  the  pressure  distributions  in 
Figs.  18  and  19  make  it  apparent  that  this  is  approximately  the  case. 

In  summary,  it  may  be  said  that  the  flow  observations  suL>.stan- 
tiate  the  analytical  models  in  that  all  three  (mixing,  sheet  flow,  expansion) 
types  of  flow  are  present.  It  is  not  possible  to  obtain  a  quantitative 
measure  of  the  extent  of  each  from  the  schlierens  but  relative  density 
gradients’^  indicate  that,  qualitatively  speaking,  the  sheet  flow  is  the 
dominant  type.  This  seems  to  agree  with  the  general  trend  of  the 
throttling  data. 

F.  CORRELATION  OF  NSL  EXPERIMENTAL  RESULTS 

In  this  section  the  Phase  II  data  is  compared  to  the  results  given 
by  the  analytical  flow  models.  This  correlation  is  discussed  further 
in  light  of  the  Phase  II  flow  observations. 

It  can  be  seen,  for  e.xample,  that  in  Figs.cTd  and  24e  that  the  sheet  flow 
portion  is  black,  rather  than  one  black  and  one  white  on  both  the  upper 
and  lower  injection  flows  in  the  Schlieren.  This  only  occurs  when  the 
density  gradient  is  so  strong  that  light  deflection  becomes  excessive 
(white  areas  become  black). 
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Some  Phase  I  results  are  compared  to  those  of  Plase  II  and  both 
are  compared  toavailable  data  from  other  sources. 

1.  Phase  II  Data  Compared  to  Analytical  Results 

The  flow  throttling  of  Phase  II  configurations  N,  Q,  R,  S,  and 
U  was  calculated  for  the  secondary  mixing,  sheet  flow,  and  secondary 
expansion  theoretical  models.  In  the  calculations,  as  explained  on 
page  13,  the  ratio  of  injection  throat  area  to  nozzle  throat  area’!'  is 
multiplied  by  an  orifice  coefficient  C  so  that  the  theoretical  and  experi¬ 
mental  stagnation  pressure  ratios  (Pog/Po)  approximately  the  same. 

The  experimental  data  and  analytical  results  are  shown  in  Figs.  27 
through  31.  It  is  evident  that  the  data  falls  primarily  between  the  sheet 
flow  and  secondary  mixing  curves.  It  was  noted  in  the  preceding  section, 
however,  that  all  three  types  of  flow  are  present  according  to  schlieren 
observations.  Nevertheless,  an  explanation  cannot  be  given  for  how 
the  flows  are  quantitatively  proportioned  to  yield  the  experimental 
throttling  results.  It  is  very  probable  that  the  proportions  are  some¬ 
what  different  for  two-  and  three-dimensional  flows.  Tliis  is  so  because, 
in  contrast  to  the  two-dimensional  data,  the  three-dimensional  results 
fall  between  the  theoretical  curves  for  the  sheet  flow  and  secondary  ex¬ 
pansion  models.  The  situation  is  further  complicated  by  the  fact  that  the 
flow  type  proportions  may  change  with  injection  rate.  In  any  event,  it  maybe 
said  that  a  qualitative  grasp  is  within  reach  with  the  analytical  models 
giving  reasonably  good  correlation- 

Z.  Comparison  of  Phase  I  and  Phase  II  Data 

All  of  the  Phase  I  and  Phase  II  flow  throttling  data  for  an  injec¬ 
tion  angle  of  60°  has  been  represented  by  "average"  curves  in  Fig.  32. 

The  averaging  process  consisted  of  fairing  a  curve  by  eye  on  a  plot 
which  displayed  all  of  the  data.  It  is  evident  from  the  curves  that 
three-dimensional  throttling  is  consistently  more  effective  than  that 
in  the  two-dimensional  case.  However,  the  discrepancy  diminishes  with 
increasing  injection  rate  beyond  a  secondary  flow  ratio  of  0.06.  In 

In  the  Phase  II  two-dimensional  case,  this  ratio  is  2b/d  where  bis 
the  slot  height  (0.08",  see  Fig  .11)  and  d  i  s  the  nozzle  throat  height 
(1.5"). 
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spite  of  the  apparent  difference,  the  general  chcaracter  of  flow  throttling 
is  the  same  in  two-  and  three-dimensional  flow, 

_ The  obvious  factor  which  may ~at  least  partially  explain  the  differ¬ 
ence  in  flow  throttling  behavior  is  the  basic  difference  between  two- 
and  three-dimensional  flow.  That  is,  the  injection  arrangement  in  the 
three-dimensional  case  may  yield  a  better  distribution  of  secondary  flow 
in  which  losses  are  smaller.  Thus,  the  secondary  flow  may  be  more 
efficiently  utilized  for  absorbing  primary  stream  momentum  and  throttling 
the  primary  flow. 

Another  factor  which  may  contribute  to  the  difference  was 
mentioned  before  and  is  the  proportion  of  secondary  mixing,  sheet  flow, 
and  secondary  expansion.  It  seems  reasonable  that  different  proportions 
might  be  associated  with  two-  and  three-dimensional  flows. 

A  small  part  of  the  difference  in  flow  throttling  might  be  due  to 
the  effect  of  injection-to-nozzle  throat  area  ratio  (A,..^/A,jj).  In  the  two- 
dimensional  case  this  ratio  was  about  twice  that  in  the  three-dimensional 
one-  The  results  in  Ref.  9  showed  that  flow  throttling  tends  to  be  more 
effective  at  lower  injection-to-nozzle  throat  area  ratio.  However,  this 
effect  is  not  pronounced  at  high  overall  nozzle  pressure  ratios  and  in 
the  present  instance  could  be  only  a  minor  part  of  the  difference  between 
two-  and  three-dimensional  throttling. 

3.  Comparison  of  NSL  and  Other  Data 

Comparative  flow  throttling  data  obtained  by  Rocketdyne  (Ref-  5) 
and  the  NASA  (Rrf.  10)  are  shown  with  some  NSL  Phase  I  data  in  Fig.  33. 
The  injection  angle  of  90"  is  common  to  all  three  configurations  although 
the  Rocketdyne  data  corresponds  to  asymmetrical  injection.  The 
Rocketdyne  data  would  have  shown  very  slightly  more  effective  flow 
throttling  had  injection  been  symmetrical.  In  view  of  experimental 
uncertainties  it  is  evident  that  the  agreement  between  the  various  data 
is  good.  Ail  of  the  configurations  in  question  are  three-dimensional. 
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SECTION  III 

THRUST  VECTOR  CONTROL 


The  majority  of  studies  of  secondary  injection  have  been  concerned 
with  the  use  of  this  method  as  a  means  of  obtaining  thrust  vector  control. 

In  the  Phase  I  work  done  by  the  NSL,  asymmetric  throat  injection  was 
investigated  to  determine  vector  control  characteristics.  It  was  found  that 
the  side  force  (component  of  the  thrust  vector  normal  to  the  nozale  axis) 
generated  was  about  the  same  as  the  jet  reaction  forced'.  In  the  Phase  II 
experiments  the  side  force  produced  by  asymmetric  injection  at  the  throat 
was  small  (a  fraction  of  the  secondary  jet  reaction  force)  and  opposite 
in  direction  to  the  secondary  jet  reaction  force.  These  Phase  I  and 
Phase  II  results  (comparable  to  those  of  Rocketdyne,  Ref.  5)  offer 
experimental  proof  of  the  fact  that  asymmetric  injection  for  vector 
control  should  not  be  at  the  nozzle  throat.  This  may  bei  e'xplained  quali¬ 
tatively  by  the  following  considerations.  Side  force  will  be  maximized 
by  tnc  production  of  a  very  asyrrmetrical  wall  pressure  distribution  with 
the  high  pressure  side  being  on  the  injection  side  (so  that  pressure  force 
adds  to  jet  reaction  force).  Much  data  has  shown  that  the  injection  effects 
do  not  extend  far  downstream  so  that  essentially  all  of  the  asymmetry  in 
pressure  distributions  must  necessarily  be  restricted  to  the  subsonic 
nozzle  walls.  However,  the  asymmetry  of  pressure  distribution  is 
precluded  to  a  large  extent  by  the  fact  that  the  flow  field  is  subsonic.  It 
is  apparent  by  this  reasoning  that  injection  at  a  supersonic  station  (with 
little  or  no  tiarottling)  will  make  the  pressure  distribution  asymmetrical 
only  in  the  supersonic  portion  of  the  nozzle.  The  degree  of  asymmetry 
will  depend  on  nozzle  geometry,  injection  rate,  and  Mach  number  at  the 
injection  station.  Injection  station  Mach  number  largely  de.termines  the 
p(5s  Hon  (if  it  exists)  at  which  the  injection-induced  shock  wave  reaches 
the  opposite  nozzle  wail.  Maximum  side  force  obviously  depends  on  the 

The  .side  force,  as  mentioned  in  Ref.  1,  was  measured  with  much 
uncertainty  so  that  only  a  qualitative  indication  was  obtained. 
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shape  of  this  asymmetrical  pressure  distribution  and  how  it  is  related 
I  to  the  nozzle  shape  (or  to  be  more  precise,  the  nozzle  area  distribution). 

Aaymmetric  injection  at  a  supersonic  nozzle  station  was  also 
investigated  in  the  Phase  II  tests  and  significant  values  of  side  force 
were  measured.  These  results  are  discussed  in  the  following  section 
while  a  correlation  with  other  data  is  presented  in  the  final  section. 

A.  NSL  SIDE  FORCE  DATA 

The  measurements  of  side  force  presented  here  correspond  to 
injection  at  the  inflection  point  of  the  low  pressure  gradient  nozzle.  Side 
force,  normalized  by  the  thrust  for  zero  injection,  is  plotted  in  Fig.  34 
as  a  function  of  secondary  flow  ratio.  The  nozzle  area  ratio  (A^jjj/A^)  at 
the  injection  station  is  1 . 85  which  corresponds  to  a  local  Mach  number 
of  2.  1  according  to  one-dimensional  flow  equations.  The  side  force 
ratio’>  is  a  linear  function  of  secondary  flow  rate  at  low  flow  rates. 
However,  it  then  levels  off  and  decreases  toward  zero  as  secondary  flow 
is  increased.  This  behavior  is  due  to  the  impingement  of  the  injection- 
induced  shock  wave  on  the  wall  opposite  injection.  The  change  in  flow 
with  injection  angle  is  seen  to  be  significant  since  more  side  force  is 
generated  at  the  lower  injection  angle  configuration  in  spite  of  a  smaller 
jet  reaction  force.  The  jet  reaction  forces  shown  are  only  the  components 
in  the  side  force  direction  (normal  to  the  nozzle  axis).  The  values  of 
aide  force  ratio  at  the  high  secondary  flow  rates  are  of  little  significance 
since  shock  "interference"  on  the  wall  opposite  injection  has  nullified 
the  effectiveness  of  the  configuration.  Injection  ac  even  higher  flow 
rates  v.'ill  not  improve  matters  since  the  induced  shock  wave  •udll  become 
more  normal  (with  nozzle  choking  po.s.sible,  eventually)  and  an  essentially 
symmetrical  pressure  distribution  which  produces  little  side  force  will 
be  maintained. 

A  series  of  schlieren  photographs  of  asymmetric  injection  flow 
patterns  is  shown  in  Fig.  35.  The  flow  configuration  near  the  wall  with 
injection  is  peneraJ.ly  similar  to  that  in  the  case  of  symmetric  injection 
(Fig.  26).  The  process  whereby  side  force  is  reduced  with  increasing 
secondary  flow  is  quite  apparent  in  Fig.  35.  The  injection-induced  shock 

The  uncertainty  in  Fg/Fj  is  ±0.005. 
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wave  impinges  on  the  wall  opposite  th6  injection  side  of  the  noazle-  The 
higher  pressure  generated  opposes  the  high  pressure  on  the  injection 
side  and  a  force  is  produced  in  opposition  to  the  side  force-  As  secondary 
flow  is  increased  this  effect  attenuates  the  side  force  further  as  the 
injection-induced  shock  moves  farther  upstream.  Losses  in  side  force 
are  also  due  to  the  decrease  in  pressure  downstream  of  the  injection 
station  at  large  secondary  flow  rates. 

An  amplification  factor  defined  (in  Ref.  5)  as  the  side  force 
ratio  (Fg/Fo)  divided  by  the  secondary  flow  ratio  (wg/wp^)  was 
calculated  for  the  a'synmctric  injection  configurations.  A  typical  plot 
of  this  amplification  factor  as  a  function  of  secondary  flow  is  shown  in 
Fig.  36.  Tlus  factor  is  a  maximum  at  low  injection  rates  but  falls  off 
rapidly  with  increasing  injection  eventually  approaching  some  equilibrium 
value  relatively  near  zero.  It  is  also  a  weak  function  of  injection  angle 
with  performance  being  slightly  improved  with  decreasing  injection 
angle.  The  variation  of  the  amplification  factor  is  obviously  a  function 
of  nozzle  shape,  exit-to-throat  area  ratio,  and  injection  station.  Best 
performance  will  be  with  a  configuration  where  the  fall-off  of  K  is  slow 
and  takes  place  at  a  high  injection  rate. 

Fig.  37  shows  axial  thrust, normalized  by  the  thrust  without 
secondary  injection,plotted  against  secondary  flow  ratio.  The  thrust 
increases  in  direct  proportion  to  secondary  flow  because  a  portion  of 
the  momentum  of  the  secondary  flow  is  in  the  axial  direction.  That  is, 
the  ratio  of  secondary  stagnation  pressure  to  nozzle  exit  pressure  de¬ 
termines  the  ultimate  velocity  and  some  of  the  corresponding  change 
in  momentum  appears  as  pressure  on  the  nozzle  wall  and  thus  axial 
thrust.  A  portion  of  this  momentum  is  negative  thrust  due  to  the 
direction  of  the  injection  slot.  The  increase  in  thrust,  as  shown  in 
Fig.  37,  is  an  increasing  function  of  injection  angle  as  would  be 
expected. 

B.  CORRELATION  OF  NSL  AND  OTHER  DATA 

The  side  force  measurements  discussed  in  the  preceding  section 
are  shown  in  Figs.  38  and  39  along  with  comparable  data  from  other 
investigations.  The  data  corresponds  to  references  as  follows: 
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UAC  {Rei.  6),  NOTS  (Ref.  11),  Rocketdyn©  (Ref.  5),  and  NASA  (Ref.  12). 

The  nozzle  exit-to-throat  area  ratios  {A^/A^),  injection  angles  (4)),  and 
ilij^ctibh  statioiis  (Ainj/A,,,j  are  seen  to-vaxy  over  a  wide  range  so  that 
a  high  degree  of  quantitative  correlation  cannot  be  expected.  The 
qualitative  agreement,  however,  is  good  for  several  reasons.  The 
NSL  and  NOTS  data  agree  at  low  secondary  flow  rates  and  examination 
shows  that  they  are  closest  in  Ag/A^,  and.  Ajnj/A,^.  No  comparison  is 
possible  at  higher  injection  rates  because  of  the  lack  of  NOTS  data.  This 
is  also  partially  true  fOr  the  Rocketdyne  and  NASA  measurements. 

However,  the  Rocketdyne  and  NASA  data  does  not  show  the  loss  in  effect¬ 
iveness  that  is  apparent  in  the  NSL  data.  Nevertheless,  the  Rocketdyne 
data  at  least  shows  a  small  decrease  in  slope  with  increasing  secondary 
flow.  It  is  believed  that  in  the  Rocketdyne  and  NASA  cases  (and  in  the 
NOTS  one  also)  serious  pressure  forces  have  not  begun  to  develop  due 
to  injection-induced  shock  impingement  on  the  wall  opposite  injection. 

This  is  probably  due  to  the  use  of  bell- shaped  rather  than  conical 
(or  slowly  diverging)  nozzles  which  tend  to  be  short  and  thus  less  vulner¬ 
able  to  shock  capture  by  the  opposite  wall. 

In  the  UAC  experiments,  however,  the  data  is  obviously  similar 
in  showing  the  onset  of  side  force  nullification.  Schlieren  photographs  are 
also  shown  (Ref.  6)  from  these  tests  which  indicate  the  actual  location  of  the 
shock  impingement.  In  spite  of  the  large  Ae/A,^,  shock  impingement 
took  place  because  the  nozzle  was  long  owing  to  a  conical  shape  and  15' 
semi-vertex  angle.  The  smaller  slopes  (at  low  secondary  flow  rates) 

.if  the  UAC  curves  are  probably  due  to  a  less  effective  injection  angle 
(less  in  terms  of  induced  flow  configuration  only^ since  the  side  force 
component  of  jet  reaction  force  is  the  same  for  UAC  and  NSL  cases 
because  sin  75°  =  sin  105°)  and  the  difference  between  two-  (NSL)  and 
three-dimensional  (UAC)  flow. 

A  graphical  tnethod  for  assessing  the  relative  merits  of  vector 
control  schemes  has  been  reported  in  Ref.  13.  This  method  compares 
the  side  force  ratio  (Fg/Fo)  to  a  specific  impulse  decrement  parameter 
(Al/I)  as  shown  in  Fig.  40.  On  the  assumption  that  the  injectant  is  equi¬ 
valent  to  the  main  (primary)  flow,  the  loss  in  specific  impulse  is 
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where  the  increment  in  thrust  due  to  injection  is  of  the  order  of  Fg  sin  i{< 
and  is  the  noiizle  wall  angle’*  (at  injection  station)  relative  to  the  axis. 
Since  Ref.  13  shows  that  range  losses  are  directly  proportional  to  loss 
in  specific  impulse,  the  most  desirable  means  of  control  is  one  with  a 
maximum  slope  which  extends  to  a  large  ordinate.  The  shapes  and 
locations  of  the  curves  in  Fig.  40  are  about  the  same  (relative  to  one 
another)  as  in  Fig.  38  and  39*  The  Hausmann  data  from  UAC  is  viewed 
with  caution  as  the  shroud  involved  may  have  been  asymmetrical  which 
could  markedly  influence  the  data  (more  favorable  side  force  with 
more  nozzle  area  on  the  high  pressure,  injection  side).  It  should  be 
noted  that  although  the  hot  firing  data  shows  poor  performance  it 
corresponds  to  a  very  high  A^^j/A^^.  The  performance  of  the  config¬ 
urations  in  Fig.  40  should  be  regarded  merely  as  a  range  of  possibilities. 
That  is,  they  should  be  viewed  as  establishing  the  characteristics  of 
the  method  and  the  influence  of  the  important  variables.  The  maximum 
potential  of  the  injection  method  should  be  evaluated  on  the  basis  of 
performance  of  an  optimization  study  design.  Furthermore,  in  the 
final  analysis  all  schemes  (injection  and  otherwise)  must  also  be  judged 
in  view  of  their  relative  weight  penalties  and  control  complexity. 


It  is  believed  that  tan  ip  should  be  used  rather  than  sin  ip  when  Fg  is 
measured  in  the  direction  normal  to  the  axis.  Ref.  13  quotes  the  sin  ip 
as  shown  above  but  for  the  range  of  vp  of  interest  here  the  difference 
between  sin  vp  and  tan  ip  is  not  important. 
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SECTION  IV 
CONCLUSIONS 

A  study  of  rocket  thrust  control  by  gas  injection*  has  been  com¬ 
pleted.  The  conclusive  results  of  this  investigation  form  the  basis  of 
the  remarks  below.  Tha  material  in  this  section  breaks  down,  as  it 
did  in  the  investigation  itself,  into  two  parts:  thrust  modulation  and 
vector  control. 

Thrust  Modulation 

The  following  items  characterize  the  method  of  gas  injection  as 
a  means  of  thrust  modulation. 

1 .  Thrust  varies  directly  with  total  nozzle  flow.  In  many  cases 
the  proportionality  is  the  same  as  that  in  the  nozzle  with  no  injection. 

In  some  two-dimensional  cases  the  proportionality  factor  increases 
slowly  with  injection  rate  while  certain  three-dimensional  ones  show 
a  similar  decrease.  These  changes  from  the  zero  injection  proportion¬ 
ality  approach  (linearly)  values  somewhat  under  10  percent  at  an  injection 
rate  of  ws/wp^  s  0.3.  The  proportionality  is  further  a  function  (weak) 
of  only  nozzle  geometry  with  sensitivity  being  atteunated  by  increasing 
expansion  ratio. 

2.  An  analysis  of  flow  throttling**  was  devi sed  based  on  three 
models .  Schlieren  photographs  verify  the  fact  that  all  three  types 

of  flow  Occur,  The  secondary  mixing  and  sheet  flow  models  correlate 
tlu-  two-dimensional  flow  throttling  data  whi’e  the  sheet  flow  and  secondary 
expansion  models  correlate  the  tiiree-dimensional  flow  throttling  data, 

Tlic  reduction  in  primary  flow  in  the  three-dimensional  case  is  about. 

2  j  times  the  injected  secondary  flow  while  in  the  two-dimensional 
case  the  reduction  is  somewhat  less.  The  constant-area  duct  assumption 
in  the  analysis  was  verified  by  schlieren  photographs  and  pressure 
distribution  data  wiiich  showed  that  the  throttling  process  takes  place 

Uriless  otherv/is'j  .snecified,  injection  is  at  the  nozrAe  throat. 

Flow  throttling  means  reduction  in  primary  flow. 
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in  a  short  axial  distance.  In  addition,  experimental  data  supports  the 
analysis  indication  that  flow  throttling  is  maximized  by  injection  at  the 
throat.  Since  thrust  varfes  with  total  nozzle  flow,  a  knowledge  of  the 
flow  throttling  implies  a  knowledge  of  the  thrust. 

3.  Flow  throttling  is  insensitive  to  orifice  geometry.  That  is, 
a  given  secondary  injection  orifice  area  may  be  distributed  as  a  slot 
or  holes  of  various  sizes  and  flow  throttling  is  unchanged.  It  is 
believed  that  shaped  injection  ports  (supersonic  nozzles)  will  not  signifi¬ 
cantly  improve  flow  throttling;  the  injection  stagnation  pressures  required 
for  filling  an  injection  port  nozzle  of  large  Mach  number  (M  >  2)  fast 
become  too  high  to  be  practical. 

4.  Primary  nozzle  geometry  has  a  small  effect  on  flow  throttling 
which  ^  attenuated  by  increasing  overall  nozzle  pressure  ratio.  Flow 
throttling  is  increased  (at  a  given  w^/wp^ )  as  the  pressure  gradient 
(implies  throat  shape)  is  decreased.  The  nozzle  shapes  (i.e.,  bell¬ 
shaped)  presently  in  use  have  a  high  8p/8x  in  the  throat  region  and  thus 
are  unfavorable  in  this  respect  without  modification  (which  would  be 
major  and  perhaps  not  worth  the  small  gain). 

5 .  Primary  pressure  ratio  does  not  affect  flow  throttling  ^ 
pressure  ratios  greater  than  about  _3.  This  means  that  for  probably  all 
rocket  engines  (exit  Mach  numbers  greater  than  about  1.5)  primary 
pressure  ratio  is  unimportant  as  regards  gas  injection  flow  throttling. 

6 .  Flow  throttling  var^e^  with  injection  angle  (6)  primarily  due 
t^  axial  momentum  contribution  of  the  injected  £a£.  The  flow 

throttling  analysis  (with  which  experimental  data  seems  to  agree,  at 
least  qualitatively)  indicates  that  at  a  given  Wg/wp^  the  primary  flow 
is  further  (beyond  that  for  ij)  =  90')  reduced  by  an  amount  directly 
proportional  to  the  product  of  the  secondary  flow  and  the  cosine  of 

7 .  The  circumferential  angle  (^)  over  which  injection  takes 
place  does  i^t  affect  flow  throttling  for  configurations  which  are 
symmetrical  about  a  plane  through  the  nozzle  axis . 

8 .  The  flow  throttling  analysis  indicates  that  performance  i s 
improved  ^  using  ^  injectant  which  possesses  a  low  specific  heat  ratio 
(y)  and  a  low  molecular  weight. 
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The  increase  in  performance  conies  from  the  increase  in  the  axial 
momentum  component  (which  opposes  primary  stream  momentum) 
of  the  injected  gas.  In  the  same  manner,  flow  throttling  is  increased 
^  increasing  the  stagnation  temperature  of  the  secondary  gas. 

9.  All  of  the  theoretical  and  experimental  study  is  baaed  on 
the  maintenance  of  &  constant  primary  stagnation  pressure  regardless 
of  Uie  amount  of  flow  throttling.  This  implies  that  primary  supply 
throttling  takes  place  upstream  of  the  nozzle  stagnation  chamber  as  it 
indeed  did  in  the  experiments.  However,  such  a  capability  in  the 
rocket  engine  is  hard  to  imagine  in  a  solid  propellant  one  and  not 
presently  attainable  in  liquid  propollant  ones.  The  turbopump  compressors 
in  liquid  propellant  engines  pump  at  a  constant  flow  rate  to  a  constant 
total  head  with  performance  being  variable  only  over  a  small  range.  There¬ 
fore,  the  use  of  gas  injection  ^  a  means  of  flow  throttling  is  ^  reality 
contingent  upon  a  flow  throttling  capability  upstream  of  the  combustion 
chamber  which  can  maintain  a  constant  chamber  pressure. 


Vector  Control 

The  following  items  characterize  the  method  of  gas  injection  as 
a  means  of  vector  control. 

1 .  The  asymmetric  injection  position  for  maximum  side  force 
is  in  Uie  supersonic  portion  of  the  nozzle.  Asymmetric  injection  at 
the  throat  is  inefficient  since  the  injection-induced  asymmetry  in 
pressure  distribution  is  not  effective  in  the  subsonic  region  and  falls 
off  rapidly  with  distance  in  the  supersonic  region  so  that  the  area 
affected  is  small  (near  the  throat).  As  injection  station  is  moved  toward 
the  throat  from  the  exit,  side  force  increases  until  enough  of  the  injection- 
induced  shock  is  "captured"  by  the  opposite  wall  to  produce  a  signirU:ant 
counter-force  which  opposes  the  side  force.  The  optimum  position  will 
thus  be  somewhere  between  these  extremes  and  a  strong  function  of 
jiozzle  geometry-  The  short,  bell-shaped  nozzles  in  use  in  present 
rocket  engines  might  be  favorable  since  shock  ."capture"  may  be  attenuated. 
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2.  Side  force  varies  in  direct  proportion  ^  Becondary  flow  rate 
before  the  onset  of  serious  shock  capture  (up  to  wg/wp^  of  about  0.  06  - 
in  configurations  tested  to  date). 

3.  The  ratio  of  aide  force  ratio  (Fg/Fo)  to  secondary  flow  ratio 
(ws/wpg)  is  maximum  ^  low  injection  rates  and  numerically  about  2_. 

4 .  Side  force  varies  with  injection  angle  (i)  through  the  injection- 
induced  pressure  distribution  and  the  reaction  force  component  (which  part 
varies  as  sin  normal  to  the  nozzle  axis.  Since  the  induced  pressure 
distribution  becomes  more  favorable  with  decreasing  ^  but  the  reaction 
component  falls  off  as  sin  ^  the  optimum  injection  angle  will  be  less 

than  90 

5.  Injection  porta  shaped  for  supersonic  flow  increase  the 
effectiveness  of  a  configurati on  m  producing  side  force. 

6 .  ^  is  believed  that  distribution  M  injection  area  over  a 
circumferential  angle  (0)  will  net  only  very  small  gains ,  ij[  any. 

7.  The  optimum  gas  injection  configuration  ^s  regards  ij, 

and  injection  station  will  depend  on  the  particular  geometry  of  a  nozzle. 

It  ^  believed  that  increasing  injected  gas  momentum ,  ^  the  use  of  ^ 
high  temperature  gas  with  a  low  x  low  molecular  weight,  will 
increase  the  resultant  side  force  (at  a  given  Wg/wp^ ). 
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'igure  1.  Flow  throttling  results  for  the  NSL  and  Modified 
Martin  models 
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SOLID  CURVES  INDICATE  SECONDARY  EXPANSION  MODEL 
DASHED  CURVES  INDICATE  SHEET  FLOW  MODEL 


*  .7 

■v 


O 

a: 


J  ).0 


D 

CONFIGURATION  F 
Aii,/wR|*  »  .060 
C  t  .74 

■-  60° _ 


PRIMARY  GAS  =  AIR 
secondary  gas  =  AIR 


CONFIGURATION  J 
A»s/i7R®  =  .105 
C  =  .69 
6  -  60® 


.02  .04  .06  .08  10  .12  .14 

secondary  flow  RATIO,  w./w, 

»  *^0 

Figure  Z.  Typical  Phase  I  flow  throttling  data 
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Figure  3-  Schlieren  photograph  of  iajectioa  into  a  supersomc  stream- 
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Figure  14-  Flow  throttliDg  of  tlie  low  pressure  gra-dieut  nozzle  at  various 
stagnation  pressures  =  60®) 


CONFIDENTIAL 


TR  448 


66 


CONFIDENTIAL 


SYM  Pq  CONFIG 


CONFIDENTIAL 


1 

i 


TR  448 


68 


CONFIDENTIAL 


CONFIDENTIAL 


TR  448 


69 


CONFIDENTIAL 


CONFIDENTIAL 


(eiM)  MnStSNd 


TR  448 


70 


CONFIDENTIAL 


Figure  20a.  Effect  of  pressure  gradient  on  flow  throttling  at  an  injection  angle  of  60° 


Figure  20b-  Effect  of  pressure  gradient  on  flow  throttling  at  an  injection  angle  of  60 
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Figure  21b.  Effect  of  pressure  gradient  on  flow  throttling  at  an  injection  angle  of  30 
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'igure  24a.  Schlieren  photograph  of  injection  in  the  high  pressure 
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Fig-.ire  24b.  Schlieren  photograph  of  injectioii  in  the  high  pressure 
gradient  nozzle. 
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Figure  24e-  Schlieren  photograph  of  injection  in  the  high  pressure 
gradient  nozzle- 


Figure  25.  Schlieren  photograph  on  injection  in  the  low  pressure 
gradient  nozzle. 
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Figure  26a.  Schlieren  photograph,  of  injection  at  a  supersonic  station 
in  the  low  pressure  gradient  nozzle. 


Figure  26c.  Schlieren  photograph  of  injection  at  a  supersonic  station 
in  the  low  pressure  gradient  nozzle. 
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'igure  26e-  Schlieren  photograph  of  injection  at  a  supersonic  station 
in  the  low  pressure  gradient  nozzle. 
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Figure  35b.  Schlieren  photograph  of  asymmetric,  supersonic-fitaUon 
injection  in  the  low  pressure  gradient  nozzle- 
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APPEWDIX  A 

METHOD  OF  SOLVING  FLOW  MODEL  EQUATIONS 

The  final  equations  for  all  of  the  flow  models  were  solved  by  a 
graphical  technique.  As  shown  in  Sections  Al ,  A2  and  A3,  a  problem 
solution  finally  derives  from  a  pair  of  simultaneous  equations  with 
unknown  (Rj/R,)*  and  M, .  However,  these  equations  cannot  be  solved 
explicitly  in  the  form  given.  Thus,  a  graphical  method  was  used  in 
which  each  of  the  simultaneous  equations  was  solved  for  (R^/Rj)*  with 
both  results  being  plotted  as  a  function  of  Mj,  as  shown  in  the  sketch 
below . 


PRIMARY  FLOW  MACH  NUMBER,  M, 


The  intersection  point  is  a  solution  with  the  corresponding  Mj  and 
(Rj/Ri)^  satisfying  the  simultaneous  equations.  A  solution  was  thus 
determined  for  each  assumed  secondary  flow  rate- 

The  area  ratio  (R^/R,)*  found  from  the  above  procedure  is  a 
flow  throttling  parameter  in  all  of  the  models.  The  continuity  equation 
(Eq.  (5))  for  the  secondary  mixing  model  shows  that 
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Thus,  the  primary  flow  throttling,  the  ratio  Wp/wp^  ,  is  determined  by 
subtracting  the  associated  value  of  Wg/v'p^  from  (Rz/Rj)*  •  In  the  sheet 
flow  and  secondary  expansion  models  £q.  (10)  shows  that 


since  the  primary  and  secondary  streams  flow  without  mixing. 
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APPENDIX  B 

SHEET  FLOW  MODEL  DETAILS 


The  isentropic  flow  area  ratio  is  a  £or;m  o£  the  continuity 
equation  for  a  perfect  gas  and,  as  shown  in  Ref.  3,  is  given  by 


A  _  1 


2 

yn 


V  +  1 

2(7^ 


where  M  is  the  Mach  number  at  the  station  where  A  is  the  area,  in  a 
streamtube  flow. 

The  value  of  M^^  in  Eq.  (10)  is  derived  from  the  assumption 
that  Pig  =  Pa  •  It  is  obvious  that 

^ 

Po  Po 


A  value  of  po  /po  is  implied  by  the  assumption  of  a  Wg/wp.  and  thus 

Zi  ffs  -  fis 

P',  Po  ’  POg 

which  implies  the  value  of  M^^  according  to  Eq.  (3). 


TR  448 


109 


CONFIDENTIAL 


vtiiaessissiii 


APPENDIX  C 


SECONDARY  EXPANSION  FLOW  MODEL  DETAILS 
The  area  ratio  Aj  /A.  ,  which  determines  Mi  ,  is  determined 

*  8 '  V  s  *  8 


Azg  n  (R,'  -  R,*) 

- ^ - 


Since  Aj,  /A.  ,  and  tlius  M^  ,  depends  on  (R2/Ri)^  ,  a  trial  value  is 
needed  for  the  (Rj/Ri)*  in  the  above  equation  in  order  to  solve  Eq.  (12) 
for  the  same  parameter.  However,  as  mentioned  on  page  8,  the  area 
ratio  (Rj/Ri)'^  found  from  the  continuity  and  momentum  equations 
(Eos.  (8)  and  (12))  is  determined  for  various  assumed  values  of  Mj  . 
Thus,  in  the  present  case  the  assumed  Mj  yields  an  (R2/R,)^  from 
Eq.  (8)  which  is  then  used  as  a  trial  value  in  the  above  equation  .so 
that  Eq.  (12)  may  also  be  solved  for  (R2/Rj)*  .  The  M2^  corresponding 
to  the  A2  /A,  in  the  above  equation  is  used  to  obtain  the  presoure 
ratio  P2g/P0g  in  Eq.  (12). 
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